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THE d STEEPEST DESCENT METHOD AND THE ASYMPTOTIC BEHAVIOR OF 

POLYNOMIALS ORTHOGONAL ON THE UNIT CIRCLE WITH FIXED AND 

EXPONENTIALLY VARYING NONANALYTIC WEIGHTS 

K. T.-R. MCLAUGHLIN AND P. D. MILLER 



Abstract. We develop a new asymptotic method for the analysis of matrix Riemann-Hilbert problems. 
Our method is a generalization of the steepest descent method first proposed by Deift and Zhou; however 
r^ ' our method systematically handles jump matrices that need not be analytic. The essential technique is to 

^ ' introduce nonanalytic extensions of certain functions appearing in the jump matrix, and to therefore convert 

the Riemann-Hilbert problem into a d problem. We use our method to study several asymptotic problems of 
polynomials orthogonal with respect to a measure given on the unit circle, obtaining new detailed uniform 
convergence results, and for some classes of nonanalytic weights, complete information about the asymptotic 
behavior of the individual zeros. 



1. Introduction 



< 

u 

C^ . 1.1. Asymptotic analysis of Riemann-Hilbert problems. The steepest descent method for asymptotic 

analysis of matrix Riemann-Hilbert problems was introduced by Deift and Zhou in 1993 |14) . A matrix 
Riemann-Hilbert problem is specified by giving a triple (S, VjA/") consisting of an oriented contour E in the 
complex 2-plane, a matrix function v : S ^ SL{N) which is usually taken to be continuous except at self- 
intersection points of S where a certain compatibility condition is required, and a normalization condition TV 
^b ' as z — > oo. If S is not bounded, certain asymptotic conditions are required of v in order to have compatibility 

00 . with the normalization condition. Consider an analytic function M : C \ S ^ SL{N) taking continuous 

■^ ' boundary values M+(z) (respectively M_(z)) on E from the left (respectively right). The Riemann-Hilbert 

^§ , problem (E, v,A/') is then to find such a matrix M(z) satisfying the normalization condition A/" as z — > cx) 

and the jump condition M+(z) = M„(z)v(z) whenever z is a non-self-intersection point of S (so the left and 
f~^ . right boundary values are indeed well-defined) . The steepest descent method of Deift and Zhou applies to 

certain Riemann-Hilbert problems where the jump matrix v(z) depends on an auxiliary control parameter, 
and is a method for extracting asymptotic properties of the solution M(z) (and indeed proving the existence 
2 ' and uniqueness of solutions along the way) when the control parameter tends to a singular limit of interest. 

The original method put forth in |14| bears a striking resemblance to the well-known steepest descent 
method or saddle point method for analyzing contour integrals with exponential integrands. A distinguished 
point on E is identified (analogous to a point of stationary phase) and an explicit change of variables of the 
form N(z) = M(z)t(z) where t(z) is a piecewisc analytic matrix is introduced in the vicinity of this point 
and it is observed that (i) the matrix N(z) satisfies an equivalent Riemann-Hilbert problem with a new 
contour En and a new jump matrix vn, and (ii) the jump matrix vn converges to the identity matrix in 
the singular limit of interest for all z bounded away from the stationary phase point. One therefore expects 
that a good approximation to N(z) can be constructed by an explicit local analysis near the stationary 
phase point. With the explicit local approximant N(z) constructed, one uses it in a final change of variables 
H(z) = N(z)N(z)~^ and observes (i) that H(z) satisfies a Riemann-Hilbert problem with a possibly new 
contour Eh and a new explicit jump matrix vh and (ii) that the new jump matrix vh is now uniformly 
close to the identity matrix in the limit of interest. This allows one to construct H(z) by iteration of certain 
singular integral equations that are equivalent to any given Riemann-Hilbert problem, and to show that H(z) 
is uniformly close to the identity matrix in any region bounded away from Eh- With additional work it may 
in some circumstances be shown that H(z) is close to the identity uniformly right up to the contour Eh 
(this usually requires more detailed information about the jump matrix). In this way, one obtains a formula 
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M(z) = H(z)N(z)t(z)^'^ for the solution that can be used to compute directly an asymptotic expansion of 
M(z) vahd in the singular limit of interest. 

Since the introduction of the steepest descent method for Riemann-Hilbert problems, there have been sev- 
eral key developments. In ^21 and ^^1 a technique was established in which one makes a change of variables 
involving a matrix constructed from a single unknown scalar function g{z) analytic in C\E. The transforma- 
tion modifies the jump matrix in a way involving the boundary values g±{z) taken on E. One then chooses 
relations between the boundary values oi g{z) such that the transformed Riemann-Hilbert problem becomes 
asymptotically simple. The desired conditions amount to a scalar Riemann-Hilbert problem for g{z), which 
is easily solved in many circumstances. A crucial feature of this method is that the dominant contribution 
to the solution typically comes from subintervals of the contour E of finite length rather than from isolated 
points. Here we therefore see an important difference between singular limits of matrix Riemann-Hilbert 
problems and evaluation of saddle-point integrals. In the contributing intervals the transformed jump ma- 
trix has a factorization (see (^ , (|76() , and (|183|l below) whose factors admit analytic continuation to the left 
and right of each such interval. A further change of variables based on this analytic factorization is carried 
out in lens-shaped regions surrounding each contributing interval. Ultimately, a model problem is solved 
(typically in terms of Riemann theta functions of genus related to the number of contributing intervals) and 
along with local analysis near the endpoints of the intervals a model for M(z) is built and compared with 
M(z) to obtain a Riemann-Hilbert problem for the error H(z). When the method is successful, the jump 
matrix for H(2;) is uniformly close to the identity and thus H(z) may be constructed via iteration of integral 
equations. Significantly, the conditions imposed on the boundary values of g{z) can often be viewed as the 
Euler-Lagrange conditions for a certain variational problem (see ^3] as well as ^1 ^J and Appendix ^ of 
this paper). 

A further development emerged from problems in which it was recognized that no appropriate function 
g{z) can be found relative to the given contour E. In |19j and later in |S] it was shown how analyticity of 
the jump matrix could be exploited to effectively deform arcs of the contour E to alternative locations in 
the complex plane such that the jump matrix maintains the same functional form; specific locations of the 
arcs are determined such that there exists an appropriate function g(z) as above. These selected arcs are 
the closest relatives in the noncommutative theory to the paths of steepest descent from saddle points in the 
asymptotic theory of contour integrals. The contour selection principle was also encoded into a variational 
problem in |19| . 

More recently ^| |^ IH], new techniques have been added to the framework of the steepest descent 
method that are adapted for determining the asymptotic contribution to the solution of a coalescence of 
a large number of poles in the matrix unknown (this is strictly speaking not a Riemann-Hilbert problem 
in the sense described above due to the polar singularities, however the problem is first converted into a 
standard Riemann-Hilbert problem by explicit transformations). The key idea here is to exploit certain 
analytic interpolants of given residues at the poles. 

For the fundamentally nonlinear cases in which the dominant contribution comes from subintervals of 
a contour, a central feature is that the analytical methods rely on piecewise analyticity of the given jump 
matrix and of the boundary values of the scalar function g{z) . For some of the cases of long-time asymptotics 
of integrable nonlinear partial differential equations [141 1^. as well as the recent long time asymptotic analysis 
for perturbations of the defocusing nonlinear Schrodinger equation |16| . the dominant contribution comes 
from isolated points of the contour E, and while analyticity is not fundamental, the asymptotic calculations 
proceed by an approximation argument, in which an analytic part is deformed away, and a (small) residual 
contribution is handled by technical and analytical prowess. The approximation argument is delicate and 
requires detailed analysis that depends sensitively on the geometry of the particular contour E. 

Far from being a mere pursuit of abstraction, a simple asymptotic technique that applies to Riemann- 
Hilbert problems regardless of whether the jump matrix is analytic or not would have immediate application 
in a number of important areas. For example, a unified treatment of the asymptotic theory of orthogonal 
polynomials on the real line with general nonanalytic weights would allow the resolution of universality 
conjectures from random matrix theory in the most natural and general context (see |10llllj for the analytic 
case, and [201 ^'^^ ^^ application of the d steepest descent method described in this paper to the nonanalytic 
but convex case). As another example, if it were possible to treat systematically problems with a large 
number of poles that accumulate in a very regular but nonanalytic fashion, then the important problem of 



semiclassical asymptotics for the focusing nonlinear Schrodinger equation with general nonanalytic initial 
data could begin to be addressed. 

In this paper, we present a new generalization of the steepest descent method for Ricmann-Hilbcrt prob- 
lems that applies in absence of analyticity of the jump matrix, and yet does not depend on an approximation 
argument for the jump matrix. While we believe that the ideas we will develop in this paper are useful in 
very general contexts, wc have chosen to focus on a particular application of the steepest descent method in 
order to demonstrate the technique. 

1.2. The essence of the d steepest descent method. As mentioned above, after changing variables 
using an appropriate scalar function g{z)^ the jump matrix is converted into a form that is well-suited for 
further asymptotic analysis. A common "target" form for the jump matrix in certain arcs of E is the following 
form 

/ e^'=(^)/<^ 1 \/ 1 0\/Ol\/l 

(1) v(x) = 

\ e-*«(^)/^ / \ e~*''(^)/' 1 / \ -1 / \ e*''(^)/' 1 

where x is a real parameter along the arc of E which for simplicity here we assume lies on the real axis 
(for a representation see Figure^, e > is the control parameter tending to zero in the singular limit of 
interest, and k.{x) is a strictly increasing real function of x that is related to the boundary values of g{z) 
on E. Suppose N(z) is the unknown satisfying N_|-(.t) = N_(a;)v(a;). With the assumption of analyticity 
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Figure 1 . The complex plane in the vicinity of a contour E (here coincident with the real 
axis) supporting a factorized jump matrix. 

of k{x), one may transform the Ricmann-Hilbert problem by introducing as a new unknown a matrix 0(z) 
defined in terms of N(z) by the following scheme: in some region lying on the minus side of the arc (the 
region labeled A_ in Figure^, set 

/ 1 

(2) 0{x + iy) ■.= N{x + iy) 

and in some region lying on the plus side of the arc (the region labeled A+ in Figure^, set 

/ 1 

(3) Oix + iy) ■.^N{x + iy) 

\ _giK{x+iy)/e ^ 

Elsewhere, set 0(z) = N(z). On has thus introduced two new jump contours, one on either side of the arc 
(these are the two dashed lines in Figure ^. However, the monotonicity of the real analytic function k{x) 
implies via the Cauchy-Riemann equations that the induced jump matrix relating the boundary values of 
0(z) on these two contours is exponentially close to the identity matrix in the limit e | 0. On the original 
arc, the matrix 0(z) satisfies the constant jump relation 

/ 1 

(4) 0+{x)=0^{x)\ 
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which can be subsequently dealt with in terms of special functions. In this paper, we show how this procedure 
can be carried out effectively when one does not have the option of extending k{x) from the contour S because 
it is not assumed to be an analytic function. Wc choose to extend k{x) in a way that does not assume any 
analyticity (see Q below). The price that must be paid is that the analogue of the matrix 0(2:) above is no 
longer analytic in the regions to the left and right of the arc; therefore this matrix cannot be the solution of 
any Ricmann-Hilbcrt problem. It can, however, be the solution of a matrix d problem (or more generally a 
mixed Riemann-Hilbcrt-(9 problem). It is into this framework that we extend the steepest descent method. 
This explains the terminology of the "9 steepest descent method" . 

As the fundamental contour S in this paper is the unit circle S^, we can now be very specific about 
what we mean by an extension of a nonanalytic function in this context. Suppose that f{6) is a C™~^{S^) 
function, m = 1, 2, 3, ... . Then, we define an extension operator Em ■ C""^^(S'^) -^ C(K^ \ {0}) as follows: 

(5) Emf{r,e) := |] Z^(-^ log(r))^ 

where (r, 9) are the standard polar coordinates for M^ . Note that this indeed defines a continuous extension 
to any annulus ?'_|- < r < r_ where 0<r_|_<l<r_<oo since £'m/(l,^) = f{&)- Also, since z — re'^ 
and z = re^*^, the fundamental differential operators of complex variable theory are represented in polar 
coordinates as 

oz 2 \or r oO J oz 2 \or r oO 

and therefore we see that if /(™~^)(6') is Lipschitz, then in particular it has a derivative almost everywhere 
that is uniformly bounded, and we have 

(7) dEmf{r,e) = _-^^(-^log(r))™-l 

2r (to — 1)! 

and 

(8) dEmf{r.e) = -^- ( i--Li(-*log(r))'"-i + 2 J] :L_li(_,log(r))M 

both holding for all r > and almost all 9 £ S^ (these formulae hold at every point of the plane if / is of 
class C™{S^)). It follows that dEmf{r, 9) vanishes to order to — 1 as r ^ 1 uniformly in 9. In fact, if f{9) 
is analytic for all 9, then the infinite series Eoof{r,9) converges uniformly in some annulus containing the 
unit circle r ~ 1 and represents the unique analytic extension of f{9). 

Generally speaking, Ricmann-Hilbert problems with rapidly oscillatory jump matrices arc equivalent to 
systems of singular integral equations with Cauchy kernel and rapidly oscillatory densities. Such equations 
can in principle be analyzed asymptotically \2<j\ . This approach requires delicate arguments of harmonic 
analysis. On the other hand, the d steepest descent method we will develop in this paper avoids such com- 
plicated reasoning. Indeed, by extending contour integration into integration over two-dimensional regions, 
the Cauchy kernel becomes less singular, and the analysis becomes correspondingly more straightforward. 

In the analytic case described briefly above, the asymptotic analysis is in general complicated by the fact 
that the procedure is valid in the neighborhood of certain intervals of S, and it turns out that a different 
analysis must be carried out in the vicinity of the endpoints of the intervals. The same would be expected to 
be true in the general nonanalytic case. In order to have the clearest possible presentation, we have chosen to 
describe in this paper the d steepest descent method in the context of a problem where there are nontrivial 
cases without endpoint issues, namely the asymptotic behavior of polynomials orthogonal with respect to 
weights on the unit circle. While the presence of endpoints complicates the analysis, they do not present an 
insurmountable obstruction, and the reader is referred to |2Uj for a description of the more general theory. 



1.3. Polynomials orthogonal on the unit circle. Let (f){9) be an integrablc 27r-periodic function satis- 
fying 0(6*) > for almost all 9. For two complex-valued functions f{9) and g{9) there is an associated inner 



product 

(9) (/, g)^ := ^ f f{0)W)<PiO) ^^ = ^ / /(arg(z))5(arg(z))0(arg(z))^ . 

This inner product leads to a system {pn(-z)}^o o^ orthogonal polynomials in the complex variable z: 

(10) pniz) ^ "fnz"^ + ^ c„jZ^ , whcrc 7„ > , 

and the defining relation is the orthonormality condition 

(11) (Pm,P«)0 = <5m« , Q<m,n<oo. 

Here the polynomials Pn{z) are considered as complex- valued functions of 6 by restriction to the unit circle: 
z = e*^. The constants 7„ have the interpretation of normalization constants, and the corresponding system 
{7r„(z)}^Q of monic orthogonal polynomials is defined by rescaling: 

(12) 7r„(z) = — Pn{z) , 0<n<oo. 

In 

The orthogonal polynomials satisfy recurrence relations of the form 



T^n+iiz) = Z7r„(z) + a„+iz"7r„(l/z) 

(13) ^ ^ 

z'^'^^TTn+iil/z) = a„+iZ7r„(2) + z"7r„(l/z) , 

for n — 0, 1, 2, 3, ... . Here, the complex constants {a„}^]^ are the recurrence coefficients associated with 
the weight cf) (also known as the Schur parameters or Verblunsky coefficients). By setting z = in the first 
equation of (|13|l . it is easy to see that 

(14) a„ = 7r„(0). 

For a general discussion of properties of polynomials orthogonal on the unit circle, see Chapter XH of 
Szego's monograph |25j . in which (among other things) the asymptotic behavior of 7r„(z) for n ^ oo is 
discussed. The extraction of asymptotic formulae for quantities related to the orthogonal polynomials of 
large degree with a fixed weight (f> on the unit circle is the type of asymptotic problem in the theory of 
general orthogonal polynomials (that is, beyond particular cases involving classical special functions) for 
which results have been known for the longest time. This can be traced to the fact that if ct){9)~^ happens to 
be a positive trigonometric polynomial, then there is a closed-form expression for the orthonormal polynomial 
Pn{z) that is convenient for analysis, as long as n is sufficiently large compared to the degree of 4>{9)~^. See 
|25l §11.2]. In other words, for certain special fixed weights (f>{9) the asymptotic formulae one obtains become 
exact as long as n is large enough. This leads to a general strategy for asymptotic analysis of orthogonal 
polynomials on the unit circle based on approximating an arbitrary given positive function (j>{9)~^ by positive 
trigonometric polynomials. 

The asymptotics described in the monograph of Szego are of a rather general character and hold whenever 
log{(j){0)) is an integrable real valued function. In the years since the origin of Szego's methods, there have 
been many further developments in the asymptotic theory. These developments move both in the direction 
of generalizing the class of weights for which the Szego asymptotics arc valid (perhaps in a weaker form) and 
also in the direction of trading generality of the weight for detail of the asymptotics. It seems that certain 
problems remain difficult to treat by these methods; in particular, it is difficult to verify convergence in a 
uniform sense, and it is difficult to characterize the detailed asymptotic behavior of zeros. There is a vast 
literature on this subject; we refer the interested reader to the memoir of Nevai |23| and the forthcoming 
monograph of Simon |25 . 

Beyond being a source of classical information about orthogonal polynomials on the unit circle, Simon's 
monograph describes a different viewpoint of the theory of these polynomials. Namely, Simon and his school 
have made great progress by exploiting the connection between orthogonal polynomials and spectral theory 
for operators that encode the recurrence relations that all orthogonal polynomials satisfy (see also |18|'). 
This theory is capable of establishing a number of very general results relevant to asymptotics in the limit 
of large degree. An important point is that the hypotheses required to establish results of this kind involve 
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assumptions about the asymptotic behavior of the sequence {a„}^i of recurrence coefficients. Indeed, the 
fundamental problem of spectral theory in this context is the construction of the spectral measure (j){6) dO 
from the finite difference operator involving the recurrence coefficients {o.n}'^=i- 

On the other hand, the recovery of the polynomials {Pn{z)}^^Q and of the recurrence coefficients {an}^^i 
from the spectral measure (j){9) dO is the fundamental problem of inverse spectral theory. A general approach 
to asymptotic problems in the theory of orthogonal polynomials in which the measure of orthogonality is the 
given data therefore involves the translation of the orthogonality conditions into the conditions making up a 
Riemann-Hilbert problem for sectionally analytic matrices. A Riemann-Hilbert formulation for polynomials 
orthogonal with respect to a measure on the unit circle was described in pP and follows closely the well-known 
Riemann-Hilbert formulation for polynomials orthogonal with respect to a measure on M discovered in |17| . 
In PP, and in a number of papers which followed (see, for example, I21l210)j the polynomial of degree n 
orthogonal with respect to a specific family of weights of the form (j){9) = e~"^^^' where V{6) = ^cos{9) 
was studied in the limit n — > oo. Note that this is a joint limit as the degree n of the polynomial in question 
appears in the measure of orthogonality as well; see §0]for a general discussion of such exponentially varying 
weights. For the large n asymptotics carried out in ^, and in subsequent works with this measure, as well 
as closely related measures (see, for example, [5]), analyticity of the weight (f){6) played a central role in the 
analysis. 

In [5] , Deift used polynomials orthogonal with respect to a measure on the unit circle to give an example 
of his theory of integrable operators. Specifically, he introduced a one-parameter family of positive, analytic 
functions 4>{9\t) and related solutions of Riemann-Hilbert Problem 12. II in §|5| below (with (j}{9) replaced by 
(/)(0; t)) to Toeplitz determinants. The n ^ oo asymptotic behavior of the corresponding solution M"(z; t) to 
Riemann-Hilbert Problem 12.11 then yields asymptotics for the associated Toeplitz determinants. Exploiting 
the analyticity of 0(6*; i), Deift outlined how one obtains an asymptotic description for ]V["(2:;<). The 
calculations in we will present in § |31 may be viewed as complementary to this asymptotic calculation of |S] , 
in that we will establish asymptotics for orthogonal polynomials under the much weaker assumption that 
(j){0) is a continuous function satisfying a Lipschitz condition. Furthermore, we show how the error estimates 
depend on smoothness properties of i 



1.4. Outline and summary of results. The polynomials orthogonal with respect to a weight given on the 
unit circle in the complex plane can be characterized in terms of the solution of a matrix Riemann-Hilbert 
problem in which the contour E is the unit circle. In §|2|we describe this Riemann-Hilbert problem, and 
then in § 121 and § 0] we study the singular limit in which the degree of the polynomials tends to infinity. In 
§13 we consider the weight function to be held fixed as the degree tends to infinity, while in §^we study the 
joint limit when the degree becomes large while the weight function is exponentially varied. A summary of 
the relevant logarithmic potential theory referred to in § 0] is given in Appendices ^ and El 

The key results we obtain in the fixed weights case are described in § 13.11 To the best of our knowledge, 
the uniform nature of the asymptotics we obtain is new to the field, as is our detailed characterization of 
the zeros. While there exist several classical methods available for the asymptotic analysis of orthogonal 
polynomials on the unit circle with fixed weight function 4'{9), with the degree n of the polynomial in question 
tending to infinity, this problem is the ideal context in which to introduce the d steepest descent method. 

On the other hand, the asymptotic behavior of polynomials orthogonal with respect to a varying weight 
on the unit circle, considered in § 0] is more challenging to obtain by more classical techniques. The results 
we obtain with the use of the d steepest descent method are stated in § 14.11 A point we wish to emphasize 
is that with the use of the d steepest descent method, the analysis in the varying weights case is no more 
difficult than in the fixed weights case. This fact distinguishes the d steepest descent method from more 
classical techniques. 

1.5. Notation. Throughout, we assume a fixed norm || • || on 2 x 2 matrices. For p = 0, 1, 2, . . . , we use the 
following induced norm on sufficiently smooth matrix functions F : R^ — > C2x2 

Here x and y are cartesian coordinates in R^. For p — Q, we may apply this norm on all functions F in the 
space L°°(R^), whereas for p > we may apply this norm to a subset of functions F in the space C^'^'^^iV?) 
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of functions with Lipschitz continuous mixed partial derivatives of all orders up to and including p — 1 . The 
finiteness of the norm indicates the uniformity of the Lipschitz condition. Since for functions F in the class 
CP~^'^{M.'^) the mixed partial derivatives of order p exist almost everywhere, the condition |||F|||p < oo can 
be cquivalently expressed as saying that F has all derivatives of total order at most p in the space L°°(R^). 
We also use the notation C^"^'^(M2\{o}) and Lg°(R2\{0}) for spaces of functions respectively in C^-^'Hl^^) 
and L°°(IR^) that vanish identically for | log(a:;^ + 2/^)1 large enough (that is, outside some annulus). 

For functions V{9) defined on the circle S^ (that is, V is defined for — tt < < tt), we also say that V 
is of class C''~^'^{S^) if the periodic extension of y to S M has k — 1 Lipschitz continuous derivatives, or 
cquivalently, has k derivatives in L°°(R). A suitable norm for such functions is given by 

(16) lll^ll|o,fc:= sup \V{0)\+ sup |T/W(0)|, 

-7r<0<7r -7r<e<ir 

since it is easy to establish that for all m satisfying 1 < m < k — I, 

(17) sup |y(")(0)| < (271)'=-" sup |1/('=)(0)|. 

-7r<e<7r -7r<e<7r 

If V{6) satisfies a Holder continuity condition, there exists a unique function N(z), analytic for \z\ > 1, 
decaying as z — > cx), and taking Holder continuous boundary values on \z\ = 1, such that 



(18) V{6) = iV(e'^) + Vo + N{e'0) , 

and Vb is the average value of V. Thus, A^(e**) is the negative frequency component of the Fourier series for 

v{e): 

(19) y{d)^ ^ V^e*^■^ with coefficients V, = — / V{e)e'"''^ dO , 



and we have 

°° V ■ 

(20) N{z):=y^^, for|z|>l. 

^-^ z^ 

We also introduce the function fJ : S"^ ^ R by the formula 

(21) n{e):^2'^{N{e"^)). 

Note that Vt and V are functions that are related by the Cauchy transform. 

Throughout the paper we will use a "bump" function _B : R ^ [0, 1] with the properties that B is infinitely 
differentiable, Bil) = 1 for |/| < log(2)/2, and B{1) = for \l\ > log(2). 

2. The Riemann-Hilbert Problem for Polynomials Orthogonal on the Unit Circle 
Consider the contour S illustrated in Figure |21 Let ri be a positive integer. Relative to the contour S we 




Figure 2. The contour E of the Riemann-Hilbert problem for polynomials orthogonal on 
the unit circle is the unit circle itself \z\ = 1 oriented in the counterclockwise direction. 

pose, for each n = 0, 1, 2, 3, . . . , the following Riemann-Hilbert problem for a 2 x 2 matrix M"(z). 
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Riemann-Hilbert Problem 2.1. Find a 2 x 2 matrix M"(2;) with the properties: 

Analyticity. M"(z) is analytic for \z\ ^ 1, and takes continuous boundary values M" (z), M" (z) 

as w tends to z with \z\ ~ 1 and \w\ < 1, jifl > 1. 

Jump Condition. The boundary values are connected by the relation 



(22) 



(23) 



1 



M.\{e'^) = M^_{e'' 



1 

Normalization. The matrix M"(z) is normalized at z ~ oo as follows: 

Z-" 



lim M"(z) 







Proposition 2.1. Suppose that the positive weight function (j){0) satisfies a uniform Holder condition \(j}{02) — 
0(^1)1 "^ K\62 — OiY" for some v G (0, 1] and with some K independent of 61 and 62- Then Riemann-Hilbert 
Problem \2. 1\ has a unique solution for each integer n>0, namely if n> 0, 



( 



(24) 

and if n = 0, 
(25) 



M"(2 



7r,i(z) 



-7,i 



^z"-i^„_i(l/z) 



1 

27ri 



7n-l 

' 2-Ki 



7r„ (s)s" 



-(/)(arg(s)) ds 



\ 



l(s)s- 



MO(z) = 



/l 



yo 



1 

2nxi 



E s 



(arg(s))ds 



-0(arg(s))ds 



1 



/n particular, Af{\(0) = a„ and M2\(0) = ~7n^i /o?" n > 1. Here, {7r„}^Q is t/ie sequence of monic 
orthogonal polynomials with respect to the weight 4> o-nd the inner product ((SJ, and {a„}^]^ is the sequence 
of associated recurrence coefficients (see (jl!-{|l ) while {7n}5^o *'^ ^'^'^ sequence of associated normalization 
constants. 

Proof. If ri = 0, then the Riemann-Hilbert problem is triangular with identity asymptotics and is trivially 
solved in closed form by a Cauchy integral, yielding H25|l . Thus from now on we consider n > 1. 

The uniqueness of the solution for n > 1 can be seen from the following argument. Continuity of the 
boundary values taken on E implies that the ratio of any two solutions of Riemann-Hilbert Problem 12.11 is 
an entire function of z that tends to the identity matrix as z -^ cx). Uniqueness thus follows by Liouville's 
Theorem. 

To derive (|24|l . first note that if M"(z) solves Riemann-Hilbert Problem 12.11 then the first column of 
M"(z) must be analytic throughout the z-plane. From the normalization condition H23|) it is then clear that 
Ml\{z) is a monic polynomial of degree n while M2i{z) is a polymomial of degree at most n~l (the leading 
coefficient of M2i{z) is not determined from the normalization condition alone). The jump condition for the 
second column reads 



(26) 

In particular, since 
^. Thus, 

(27) 



M^2+{en - M^2-{<in 

satisfies a Holder condition, we may express M['2(z) as a Cauchy-type integral using 



^22+{en - M^2-{^^'') 



A/{\(e''')e-'^ 
Af^\(e*^)e-™V(6l). 



MUz) 



1 

27ri 



M'^^{a)s-'^ 



3 
8 



(arg(s)) ds , 



and the normalization condition H23|l then requires that this Cauchy integral be of order z~"~^ as z -^ oo for 
each fixed n G Z_|_. Expanding the Cauchy kernel in a geometric series we see that the following conditions 
must be satisfied 



(28) i A/i"i(z)z*''^"</.(arg(z)) dz ^ Q 

for fc = 0, 1, 2, . . . , n — 1. Since dO = dz/{iz) for an angular coordinate 9 on the contour S, this proves that 
M{\(z) is orthogonal to the monomials 1, z, z^, . . . , z"~^ with respect to the inner product ©• The existence 
of such a monic polynomial of degree n follows from the Gram-Schmidt algorithm. Thus, 7\/{\(z) = 7r„(z), 
the nth monic orthogonal polynomial with respect to the weight (j){0) on the unit circle. 

A similar argument applies to the second row of M"(z). Indeed, M^j^(z) is a polynomial of degree at most 
77,-1. Using (|26|l we may express M|'2(z) as a Cauchy integral: 

(29) M^,{z) = — i ^^ 0(arg(s)) ds , 

Zm Jy: s — z 

and then the normalization condition (|23|l requires that M22{z) = z~" + 0{z~'^~^) as z — » cx). Expanding 
the Cauchy kernel in a geometric series, one sees that M2i{z) is required to satisfy the following conditions: 

(30) (Af^\(e^''),l)^ = -1 and (Af2'\(e^'^), e"='')^ = , fc=l,...,n-l. 
Equivalently, these relations may be written in the form 

(31) /e*("-i)e,e'("-i)«M2'\(e'»)\ =-1 and U^\ e''^''-^^'^lA^J^\ =0, fc = 0,...,n-2. 



Clearly, the degree n — 1 polynomial z"^ Af2\(l/z) is orthogonal to the monomials 1, z, . . . , z"^ with respect 
to the inner product (, )^, and the normalization condition then fixes the leading coefficient: z"^^7\f2i(l/z) = 
— 7^_iz"^-^ + .... In other words, we have foimd that z"^-^7\f2\(l/z) = — 7,^j_j7r„_i(z), or equivalently, 
M^^{z) = -7^iz"-i7r„_i(l/z). D 

Thus, the conditions of Riemann-Hilbcrt Problem 12.11 serve to define the orthogonal polynomials as 
an alternative to other representations that may be available, possibly including explicit contour integral 
formulae (for special families of weights). 

Riemann-Hilbert problems like this one frequently arise as a consequence of the application of the Fourier 
transform or z-transform to certain types of linear integral equations (i.e. the Wiencr-Hopf technique, see 
[7]). Reversing this sort of reasoning, the representation of the orthogonal polynomials in terms of Riemann- 
Hilbert Problem 12.11 immediately yields integral equations for certain auxiliary unknowns. Of particular 
interest are the Marchenko equations obtained in |18| . in which the unknowns are Fourier coefficients of 
functions explicitly related to the orthogonal polynomial 7r„(z), and the relevant operator is of the form 
I — K, where K is an integral operator acting in ,^^(n-|- 1, n-|- 2, . . . , oo) with a kernel that depends explicitly 
on the weight (j) but not otherwise on the degree n of the polynomial in question. Furthermore, the kernel 
docs not depend in any crucial way on the smoothness of the weight. Therefore, in principle, this formulation 
makes possible the calculation of asymptotics for the polynomials (indeed, this is one of the applications of 
the Marchenko equations discussed in [181 ) in a way that is relatively insensitive to the analyticity properties 
of the weight function (j). However, the correction terms that appear in such a scheme are necessarily in 
terms of infinite Fourier series (see, e.g., equation (VI. 7) of JHj) that while having known coefficients are 
not convenient for detailed analysis of zeros of the polynomials in regions of the complex plane where these 
zeros necessarily arise from a competition between different terms in an expansion. To provide details of the 
asymptotics, it is more advantageous to work with the Riemann-Hilbcrt problem directly. 

The uniqueness of M"(z) coupled with symmetry of the Riemann-Hilbert problem under reflection through 
the unit circle leads to the following result. 



Proposition 2.2. The matrix M."" (z) satisfying Riemann-Hilhert Prohlem \2. 1\ satisfies the symmetry relation 



(32) M"(z) = i'"^M"(0)-iM"(l/z)(-iz)'"^ 

9 



In particular, by taking z ~ above, we obtain the identity 

2 



(33) 



7«-i 

In 



1 



3. Fixed Weights 

3.1. Asymptotic behavior of orthogonal polynomials and related quantities as n — > oo. In this 
section we describe several asymptotic results in the theory of orthogonal polynomials with fixed weights 
on the unit circle that we will obtain as a fundamental illustration of our method. It will be convenient to 
introduce the real-valued function T^ : 5^ -^ R such that 

(34) (P{e) = e-"^^"^ , for aU 6 e SK 

The fundamental object of the asymptotic theory for the fixed weight (j) is the so-called Szego function: 

1 f V{arg{s))ds' 



(35) 



3^2 



exp 



27ri 



z ^S. 



This is a function analytic for |z| 7^ 1 that decays to zero as z — *■ 00. Its value at z = has the interpretation 
of the geometric mean of the weight 



(36) 



^,(0)=exp(-i-£ 



¥(9) d9] =e 



_.-Vo 



If 4'{9) satisfies a Holder continuity condition, then by strict positivity so does V{9). In this case, by the 
Plemelj formula |25], we have 



(37) 



limS'^(re*^) 

rfl 



limSJre'' 

ril 



Furthermore, recalling the negative frequency component N{z) of V defined by (|20|l . we can obtain the 
following equivalent characterization of S^{z). Since the function N{l/'z) is analytic for \z\ < 1, it follows 
that 



S^{z) 



\z\>l 
\z\<l. 



(38) 

Recalling the function VL : S^ 

(39) ^ , , 

i-Ti i-ii 

3.1.1. General theorems. The following results hold for weights where V is of class C''~^'^{S^) with k > 1. 

Theorem 3.1. Suppose that ^(6*) = e~^(^) where V : S^ ^R is of class C'''-'^'^{S^) with fc > 1. Then, for 
each fixed integer p and for each p > 1 there is a constant Kp,p > such that the estimate 



g-Vo-Nil/z) ^ 

i defined by (|21() . we have 
limS'0(re*'') • limS'0(re*'') = e-^°e'^'^^^ 



(40) 



sup 

I--I>p 



dP 



dzP 



■Kn{z)z 



-n-N(z) 



1 



< K, 



log(n) 



P.P „2fc 



holds for all n sufficiently large. 



The constant Kp^p typically blows up as p — > 1, and only a finite number of derivatives can be controlled. 
More generally, we have the following result. 



Theorem 3.2. Let p > be a fixed integer. Suppose that (l){9) = e ^'^^^ where V : S^ 
(^fe-ia^^i-j ^j^/j fc > 2p+ 1. Then there exists a constant Kp > such that the estimate 



is of class 



(41) 

holds for all n sufficiently large. 



sup 



dP 



dzP 



TTn{z)z "e 



Niz) _ ^ 



<K, 



log("-) 

P „k-2p 
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<i Remark: Note that as a special case of the estimate (|41|l we obtain the following estimate (under the 
same conditions) characterizing the polynomials on the unit circle: 



(42) 



sup 

-7r<6l<7r 



de 



^n{e'')e- 



-^Nie'" 



- 1 



<K, 



log("-) 

P „k-2p 



In fact, the proof of Theorem 13. 21 is to first establish (|42|l . from which the estimate (|41ll follows (with the 
same constant Kp) via the maximum modulus principle. > 

The weakest conditions under which the above theorem provides large-degree asymptotics are that is a 
strictly positive weight that is Lipschitz continuous. Theorem 13.21 may be compared with results reported 
in the classic monograph of Szcgo |25l §12.1]. While asymptotics of TTn{z) have been established by other 
methods under weaker conditions than Lipschitz continuity and strict positivity of the weight (/), Theorem 13. 21 
exhibits clearly the dependence of the rate of decay of the error on the smoothness of (j), and the number of 
derivatives desired. 

To our knowledge, the results of Theorem 13.21 are stronger than those previously known in that they 
establish the convergence in a uniform sense. This leads to the following. 

Corollary 3.1. Let p > be a fixed integer. Suppose that (j>{0) — e~^'^' where V : S^ ^ M. is of class 
C'=-i^i(5i) with k>2p+l. Then 



.(p)i 



(43) lim -.^II^M^I. 

n^oo nP ||7r„||0 

Proof. This follows directly from Theorem 13.21 Indeed, upon carrying out the differentiation in (|42|1 . and 
combining this estimate with its analogues for all smaller values of p, we learn that Inn {e^^)\/nP converges 
uniformly to |7r„(e*^)| for — tt < 9 < tt. The proof is then complete since on 5*^ uniform convergence implies 
convergence in L^. D 

< Remark: Notice that Thcorcm l3 . 2l immGdiatGlv implies the following formula valid for all z with \z\ < 1: 



(44) 

where 

(45) 



7r„ (z) = 



s"e^(^)+/i„(s) 



ds. 



|s| = l 



,, , ,, ^, log(n) 
sup |/i„(s)| <Rp- '^^ ' 



\s\ = l " 

While in principle this could be used to compute asymptotics for 7r„(z) in this region, more detailed analysis 
gives the following improved results. t> 

Theorem 3.3. Suppose that (/){0) = e^'^C') where V : S^ ^ R is of class C''-^'^{S^) with k>l. Then for 

> such that the estimates 

log(n) 



each p satisfying < p < 1 there are constants K^ > such that the estimates 



sup 

p<|z|<l 



(46) 
and 
(47) 
hold for all n sufficiently large. 



■^n{z) 



^«g-Vo-N(l/-)gEfcV(r,e) 



< k: 



sup |7r„(z) 

\z\<p 



K+ 



An immediate corollary is that there exists an annulus inside the unit circle that asymptotically contains 
no zeros. That the result we are about to state in this direction is in a sense sharp will be made clear when 
we consider more specific weights below in S 13.1.21 fin particular, see Corollarv l3.4|) . 

Corollary 3.2 (Zero-free regions). Suppose that (f>{9) ~ e^^^^^ where V : S^ -^ R is of class C'^~^'^(S'^) 
with k > 1. Let S > be an arbitrarily small number. Then there are no zeros o/7r„(z) in the region 



(48) 



\ogi\z\) > -{k - 5) 



log(n) 
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as long as n is sufficiently large. 



Proof. This follows immediately from the estimate (|46|) . Indeed, since iV(l/z) and Ei^V{r,d) are bounded 
for p < r < 1, zeros of 7r„(z) in the region p < \z\ < 1 necessarily arise from a balance between z" and a term 
of uniform size log(n)/n'^. However, z" is large compared with \og{n)/n'' in the region where the inequality 

log(|z|) > -{k - S) log(7i)/7i holds. D 



is of class 



A second corollary is an immediate consequence of (|47ll . 

Corollary 3.3 (Recurrence coefficients). Suppose that (j){9) — e~^(^' where V : S^ ^ 
(^fe-i.ij-^i-j ^^^^ fc > 1. Then there is a constant K > such that the bound 

(49) |a„| < ^ 

holds for sufficiently large n. 

Proof. This follows directly from H47|) with the use of the identity a„ = 7r„(0). 

Finally, we have the following result concerning the asymptotic behavior of the normalization constants. 

R is of class C''-^^^{S^) with fc > 1. Then there 



D 



Theorem 3.4. Suppose that 0(6*) = e"^^*) where V : S^ 
is a constant K > such that the bound 



(50) 

holds for sufficiently large n. 






1 < K 



log(n) 



-,2k 



We give a direct proof of this theorem based on the identity 7^_i — —M2i{0) in 5 13.3.31 However, another 
proof with a less sharp error estimate may be based upon Theorem 13 . 21 because on S^ uniform convergence 
implies convergence in L^. Thus, since ||Pn(2)||0 = 1 and Pn{z) = 7n7r„(z), 



(51) 

Using Theorem 13.21 one finds that 

(52) 

Next, using (|S5|) . we have 

(53) 



1 



rn={^l Kie^Vc^mde 



^" ^ 27r 



.Nie^") 



,N{e'^ 



dO 



^ gAf(e"')gAf(e»'') 



o 



^V(e)-Va 



log(7l) 



-Vo 



Substitution into (|52|) completes the alternate proof. 

<i Remark: At this point it is important to comment that the d method wc develop below in ii l3.2l viclds 
new formulae for the polynomial 7r,i(z) (see, for example, (|128|l ). The formulae arc semi-explicit, in that they 
arc written in terms of the solution of a 9 problem (or, cquivalcntly, in terms of the solution of an integral 
equation). This d problem is arrived at after a sequence of explicit transformations, and we prove that this 
problem has a unique solution, which possesses an asymptotic expansion for n — > oo. In general, the terms 
in this expansion can be estimated (from above). Such estimations give rise to the general results described 
in this subsection. However, in the situation that some further information about the weight function e^^ is 
known, it is frequently possible to obtain much more precise information about the terms in the asymptotic 
expansion. To illustrate what can be obtained from an analysis of the terms of the expansion, we consider in 
the following subsection a slightly more specific family of weights, and present a rather complete description 
of the pointwisc asymptotic behavior of the polynomials. > 



3.1.2. More specific weights. While the estimate (|47|l allows one to bound the recurrence coefficients, it 
does not provide an asymptotic description of the polynomial 7r„(z) for z bounded within the unit circle. 
In particular, H47|l is insufficient for deducing the location of the zeros. With further assumptions on the 
regularity of V{6) we can extract a leading term that paves the way for further analysis of 7r„(z) outside the 
zero-free region, but within the unit disk. 

Theorem 3.5. Suppose that (j){9) = e"^(^) where V : S^ ^ R is of class C'^^^^'^iS^) with k >2. Suppose 
further that V^^' (6) is piecewise continuous with ^ < oo jump discontinuities at points —tt < 6i < 62 < ■ ■ ■ < 
9i < TT, of magnitudes 

(54) aI-''^ := lim V^''^ (9) - lim V^'''> (0) . 

Let V^''' (9) have one Lipschitz continuous derivative between consecutive jump discontinuities. Then, for 
each e > 0, (7 > 0, and 5 > Q, the estimate 



(55) sup n'=+ie-^(i/^V„(z) - n''+^z''e-^»-^''^^'~^e^^^^^''^ B{\og{\z\)/e) ~ fn{z) 

log([2|)<-(fc-cr)log(ri)/n 



<<^, 



holds with 

for all n sufficiently large. 



(56) /„(,):^!_^Afe^^^(«^)^ 



Note that /„(z) is a rational function of z with poles at the £ points of discontinuity of V^''\9) on the 
unit circle, and with £ — I zeros which may lie anywhere in the complex plane, and fluctuate about as n is 
varied. 

With this result, we can completely characterize the zeros of 7r„(z) under the same assumptions on V{9). 
The simplest example of orthogonal polynomials on the unit circle is of course the case (j){9) = 1, in which 
case 7r„(z) = z" for all n > 0. Here we see that all zeros of TTn{z) lie exactly at z = 0. In particular, the zeros 
avoid the unit circle \z\ = 1. This situation is typical for strictly positive analytic weights, in which case 
it is known that the zeros of 7r„(z) asymptotically lie within a smaller disk |z| < p < 1. Here, the nearest 
singularity zq to the unit circle of the analytic continuation through \z\ = 1 of the function S^{z) from the 
domain |z| > 1 determines the radius p hy p = \zo\. However, such confinement of the zeros within the 
circle is no longer typical once one leaves the analytic class. For example, discontinuities in any derivatives 
of (j)i9) make it possible for at most a finite number of zeros to be bounded away from the unit circle while 
all remaining zeros converge to the unit circle, as the following corollaries of Theorem 13.51 show. For each 
M > let 

F+{M) :=: {zeC such that log(|/„(z)|) > M} , 
(57) F~{M) := {z G C such that log(|/„(z)|) < -M} , 

F°{M) := {z e C such that | log(|/„(z)|)| < M} . 

Corollary 3.4 (Zeros near the unit circle). Assume the same hypotheses as in Theorem \3.5l Let A„(cr) 
denote the annulus 



(58) A„(a) 



(fc + l)i^^--<log(|z|)<-(fc-.)-^°^(") 



n n 



Then, for each tr > there is some M > such that the region j4„(cr) fl F~ [M) contains no zeros of iTniz) 
for sufficiently large n. 

For each M > 0, the zeros o/7r„(z) in the region A„((t) fl {F^{M) U F^{M)) satisfy 

(59) 1^1 = 1 _ (fc + i)l2lM + ilog(|/„(z)|) + {- 

n n \n 
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and 

(60) e = --n{e) + - arg(/„(z)) + -+o(- 

n n n \n 

modulo 27r/n, where 6 = arg(z), and in both cases the error term is uniformly small in the specified region. 
It follows that the angular spacing between neighboring zeros o/7r„(z) in the specified region is A9 = 27r/n + 
o(l/n). 

For any fixed M > 0, ()59|l can be rewritten uniformly in the region v4„(a) n F^(AI) as 

(61) |zHl-(fc + l)M!)+ofi 

n \n 

and consequently there exists some a G (0, cr) such that the zeros o/ 7r„(z) in the region A„(ct) D F^{M) 
asymptotically lie between the two circles |z| = 1 — (fc + 1) log(n)/n ± a/n. 

If M > is sufficiently large, then the region F^{M) is contained in a disjoint union of small discs 
centered at the poles e ' of the rational function fn{z). In this situation, let F^ AM) denote the component 
of F^{M) near the pole e*^^ . Then from ()59|l we see that the zeros o/7r„(z) in the region v4„(cr) n F^ AM) 
satisfy 

(62) |^| = i_(fc + i)M!)_iiog|z-e'«^|+ofi 

n n \ n 

where the error term is uniform in the specified region, which indicates that zeros are attracted to a curve 
that "bulges" outward from the circle \z\ = 1 — (k + 1) \og(n)/n in a region of angular width proportional to 
log{n)/n centered at the point e*^^, to a maximum radius defined by the equation 

(63) i_|,|^fcM!)+M^lM)+ofi 

n n \n 

Note that this radius is just within the inner boundary of the zero-free annulus described by Corollary 

Proof. The annulus An{<j) converges toward the unit circle as n —>■ cx), and therefore if z is a zero of 7r„(z) 
and n is large enough, Theorem 13.51 gives 



(64) ^k+l^J^^-Vo-2Nil/z)^E,Vir,0) ^ ^^(^) ^ ^^) ^ 

as ri — > cxD because for such z, B(log{\z\)/e) — 1 and 7r„(z) = 0. The o(l) error term is uniformly small for 
all zeros in A„(cr). Now, let 



(65) c:= inf 



g-Vo-2W(l/-)gBfcV(r,e) 



and note that c > due to the assumptions in force on V . Since 
(66) inf ?i''+^\z\" ^ e-" , 

we see that ()64(l is inconsistent for large enough n if |/n(z)| < ce"'^. Therefore, given a > 0, An{(j) r\F~{M) 
contains no zeros as n — > c» as long as M > a — log(c). 

For any M > we now consider those zeros z of 7r„(z) in the region An{a-) n (F°(M) UF+(Af)), in which 
case we may divide through in H64|l by fn{z) to obtain 



(67) :^__^^-Vo^Nil/^)^E,Vir,9) ^ ^ ^ ^(^) 

fn{z) 

Consistency requires that n^^^z"^ / fn{z) = 0(1), and then since the n-independent exponential factors are 
continuous up to the unit circle and Anicr) is converging to the unit circle, we may replace these factors by 
their limiting values on the unit circle without changing the error estimate. Therefore, (|67|l becomes 

(68) !^^e^^W + 1 = 0(1), 

as n — > oo uniformly for those zeros z of 7r„(z) that lie in the annulus An{(7). This proves both (|59() and 
(EOl). D 
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o Remark: Note that the presence of the outward "bulges" in the zero curve near the points of dis- 
continuity of V^'^^ {8) indicates the sharpness of the zero-free region estabhshed for more general weights in 
Corollary El > 

While most zeros of 7r„(z) move toward the unit circle as ri -^ cx) under the hypotheses of Theorem 13.51 
there may be at most £ — 1 zeros further inside the unit circle, which correspond to zeros of fn{z). We refer 
to these as "spurious zeros" . 

Corollary 3.5 (Spurious zeros). Assume the same hypotheses as in Theorem, \S.5[ For each M > 0, there 
exists a a > 0, such that the zeros o/7r„(z) lying in the disk log(|z|) < — (fc-|- l)log{n)/n — a/n also lie in 
the set F~ (M) for sufficiently large n. Moreover, whenever e^ is a sequence oj positive numbers such that 
e„ ^ as n —^ oo, the zeros of 7r„(2:) in the disk \z\ < 1 — {k + 1) log(n)/n — l/(ne„) satisfy 

(69) /„(z) = o(l) 

as n —^ oo. In particular, iTniz) has exactly one zero for each zero of fn{z) in this region, making at most 
£ — I spurious zeros. 



sup 

\z\<l 



Proof Define a constant C > by 

(70) C -.^ 

Therefore, 

(71) sup 

log(l^l)<-(fc+l) log(n)/«-cr/n 



-2W(l/z)gBfcl/(r,( 



B{\ogi\z\)/e) 



^fc+l^„^_y„_2W(l/-)gB,y(.,e)^(lQg(|^|)/^) 



< Ce- 



lt then follows easily from H64|) that all zeros of 7r„(z) lying in the disk where the inequality log(|z|) < 
— (fe -I- 1) log(n)/n — a/n holds will also lie in the set F~{M) for large enough n whenever a > M + log(C). 
The term in H64() proportional to n'^^^z" is o(l) as n — > cxa uniformly for z in the disk delineated by the 
inequality 1^1 < 1 — (fc -|- l)log(n)/n — l/(ne„) with e„ — > as n ^ cx). Therefore, fniz) = o(l) for zeros 
in this disk, and the one-to-one correspondence of zeros of fniz) with spurious zeros of 7r„(z) in this region 
follows from the Implicit Function Theorem. D 

< Remark: Note that the zeros of fn{z) play an apparently contradictory role in the asymptotics. Indeed, 
zeros of fniz) that occur near the unit circle repel zeros of 7r„(z), while each zero of fniz) that occurs far 
enough within the unit circle attracts precisely one zero of 7r„(z). !> 

<i Remark: A careful reading of the proof of Theorem 13. 51 (see S I3.3.3|) shows that the o(l) estimate that 
is stated in H55|) can be improved to give a rate of decay, and that even better estimates can be obtained if one 
does not insist on uniformity. These simple improvements can provide, for example, decay rate information 
for the error terms in the description of the spurious zeros. We have opted not to give these slightly improved 
estimates in the interest of simplicity of presentation. i> 

3.1.3. Numerical computation of zeros o/7r„(z) when derivatives ofV have jump discontinuities. To illustrate 
the detailed asymptotic behavior of the zeros of 7r„(z) explained above, we have carried out some numerical 
experiments. Let us fix £ angles of discontinuity {6*1, . . . , 9(} C (— tt, tt) by the formula 

27r / 1 



(72) 

Consider the family of weights 



■^2 2, 
given by the formula 



for j = 1, 



J. 



(73) 



:= < 



sm - 



for 0j^i <e <9j with j = 2,3,. 



for \e\ > i£~l)TT/£. 



The positive integer k and the real numbers wi., . . . ,wi are free parameters. This weight is of the form 



e~^ where V^'^'iO) has jump discontinuities at the points 9j of magnitudes that can be adjusted by 
u-j . One advantage of this family from the point of view of numerical computation is that 
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choice of the w-j 



the Fourier coefScients of (j){0) can be evaluated symbolically. In a package such as Mathematica capable of 
arbitrary precision arithmetic, this leads to the possibility of computing the elements of the Toeplitz matrices 
(whose minors arc assembled to yield the coefficients of the polynomial 7r„(z)) with sufficient accuracy for 
subsequent numerical computation of the zeros when n is large. In practice, wc computed the coefficients 
up to an overall factor by scaling the Fourier coefficients making up the Toeplitz matrix by e^°. This is 
necessary to avoid numerical overflow or underflow since according to the strong Szcgo limit theorem the 
Toeplitz determinant of = e~^ of size n + 1 scales as e""^''. While the Fourier coefficients can be computed 
symbolically, we obtained the coefficients A^ ^e*^*^^^^ appearing in the rational function f„{z) defined by (|56|) 
with the help of numerical integration. 

The Mathematica code we wrote to carry out these computations is available from the companion website 
to this paper: "http : //www . math . Isa . mnich . edu/~millerpd/ResearchAndPublications/CP . html" The 
code takes as input the number of jump discontinuities, £, the order k of the derivative experiencing the 
discontinuites, a vector w of length £ containing the parameters Wj , and the degree n of the polynomial 7r„ (z). 
The output is a figure showing the unit circle (black) with exterior tick marks at the angles 0j, j — 1, . . . ,£, 
the zero-attracting circle |z| = 1 — (fc + l)log(n)/n (green), the inner boundary circle \z\ = 1 — k\og(n)/n 
of the zero- free annulus (red), and the zeros of fn{z) that occur in a neighborhood of the unit disk (large 
lavender dots). Superimposed on the figure arc the zeros of 7r„(z) (small black dots). Sample output from 
the program is shown in Figure |31 




Figure 3. The output of the Mathematica code for n — 104, k ~ 3, £ ~ 7, and w'^ = 
(-1, -1/2, -1/4, -1, -1/4, -1, -1/2). 

The first effect wc would like to illustrate is the rate of convergence of the zeros of 7r„(z) to the unit circle 
with increasing n. See Figure 0] The annulus associated with the inequalities 1 — fclog(ri)/ri < |z| < 1 is 
asymptotically zero-free, and the curve [zj = 1 — (/c -I- 1) log(n)/n asymptotically attracts the zeros near the 
unit circle. The convergence to the zero-attracting circle is clear. More difficult to discern from the images 
is the outward "bulging" of the zeros near the angles 0j of discontinuity toward the inner boundary of the 
zero-free region. The imaginative reader can see this effect beginning in the figure corresponding to rt = 160, 
but larger values of n (and a rcscaling of the figures near the unit circle) will be necessary to resolve the 
"bulging" completely. 

Zeros of fn{z) play little role for the polynomials whose zeros are illustrated in Figure^ Next, we would 
like to illustrate the effect zeros of fn{z) can have on 7r„(z); this is the phenomenon of spurious zeros. Note 
that in the present case of equally-spaced angles 9i, . . . ,9i, the function fn{z) is periodic in n with period £. 
In FigureElwe present images corresponding to one period of the function fn{z) in the case of discontinuities 
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Figure 4. A sequence of plots showing convergence of the zeros of nn{z) toward the curve 
|z| = 1 — (fc + 1) log{n)/n as n ^ cxd. The parameters are £ ~ 3, W"^ — (—4, —2, —3), and 
k = 2. From left to right, n ^ 10, 20, 40, 80, and 160. 

of the second and third derivatives of V{6). Here it is clear that the i ~ 1 zeros of fn{z) fluctuate about 




Figure 5. Periodic fluctuation of the spurious zeros with respect to n. Top: k ~ 2. Bottom: 
fc = 3. In both cases n varies from n = 60 to n = 66 from left to right, there are £ = 7 equally- 
spaced points of discontinuity, and W"^ ~ (—1, —1/2, —1/4, 1/2, —1/4, 1/4, —1/2). Note that 
while for k = 3 the majority of zeros appear to lie in the zero-free annulus 1 — klog(ri)/n < 
\z\ < I, this is a finite n effect. 

rapidly with n, and can be either inside the unit circle or outside. Each zero of /n(z) inside the unit circle is 
an asymptote for exactly one zero of 7r„(z), while those outside the unit circle have little effect on the zeros 
of 7r„(z). The zeros near the unit circle, either inside or outside, have a repulsive effect on the zeros of 7r„(z). 

From the images in Figure |S1 it is not obvious that the zeros of /„(z) inside the unit disk attract cor- 
responding spurious zeros of 7r„(z) in the limit n —> oo. The images shown in Figure El show that this 
convergence indeed occurs. Here, we have used the periodicity of /n(z) to examine the asymptotic behavior 
of the zeros of 7r„(z) along a periodic subsequence of ri- values along which the zeros of fniz) remain fixed. 

The effect of a zero of fniz) upon those of 7r„(z) is the most subtle when it occurs near the unit cir- 
cle. It should be stressed that the parameters Wj of the weight under consideration can be deformed 
in a continuous manner such that it may always be arranged that fniz) has zeros near the unit circle. 
We fixed n and chose a one-parameter deformation of the Wj = Wjit) in order to continuously tune a 
zero of fniz) through the unit circle from outside to inside. A movie of this deformation is available at 
[http : //www . math . Isa . umich . edu/ "millerp d/Resear chA ndPublicat ions/CP . htinl| and several consecu- 
tive frames of this movie are shown below in Figure Here it can be clearly seen that as a zero of fniz) 
enters the unit disk, it initially repels the zeros near the attracting circle of \z\ = 1 — ik -\- l)log(n)/n by 
pushing them inwards. The zeros along the attracting curve then move apart to make way for the incoming 
zero of fniz). Exactly one zero of 7r„(z) fails to get out of the way, however, and instead enters the orbit 
of the moving zero of /„(z). As the zero of fniz) moves inside the attracting circle, it thus draws with it a 
spurious zero of 7r„(z). 

3.2. The d steepest descent method for fixed weights. Here, we begin the task of proving the theorems 
stated in S I3.1l bv analyzing the behavior of the matrix M"(z) solving Riemann-Hilbert Problem 12. II for a 
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Figure 6. Convergence of spurious zeros with increasing n. Top: k = 2 and n = 61, 68, 
75, 82, 89, 96, and 103 from left to right. Bottom: /c = 3 and n = 62, 69, 76, 83, 90, 97, and 
104 from left to right. In both cases there are £ = 7 equally- spaced points of discontinuity 
and w^ = (-1, -1/2, -1/4, 1/2, -1/4, 1/4, -1/2). 




Figure 7. A continuous one-parameter deformation of a weight of the form H73() with 
e = 7 and k ^ 2. The vector of parameters is w(t)^ = (0,0,0,1/2,0,1/4,0) - 
(1, 1/2, 1/4, 0, 1/4, 0, l/2)i and t varies from t = 1/40 to t = 13/80 in steps of At = 1/80. 
The frames are ordered left to right and top to bottom. 



fixed weight 0, in the limit n ^ oo. Wc recall the representation (|34|l of 0(0) in terms of T/^ ; S*^ ^ M. In 
force, in order that Riemann-Hilbert Problem 12 . II indeed describes the orthogonal polynomials with respect 
to (p{0), is the following. 
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Assumption 3.1. V is a real continuous function on the circle that, for some exponent v e (0,1] and for 
some constant K > 0, satisfies a uniform Holder continuity condition \V{62) — V{0i)\ < K\92 ~ OiY" . 



This guarantees that (j){9) is a strictly positive function that also satisfies a Holder continuity condition 
with the same exponent, but with a possibly different constant K . 

3.2.1. Conversion to an equivalent d problem. Solution of the d problem in terms of integral equations. We 
proceed in several steps. First, let N"(2) be a new unknown related to M"(z) as follows: 



(74) 



N"(z) 







M"fz) 







for Izl > 1 



for Izl < 1 . 



[ M"(z), 

It follows from Riemann-Hilbert Problem l2.1l that the new unknown N"(z) tends to the identity matrix as 
z -^ oo, and that N"(z) is analytic for \z\ ^ 1, with boundary values on the miit circle related by 



(75) 



N"(e*'') =N!!(e*'') 







'inO 



where N"(z) (respectively N":(z)) indicates the boundary value taken at the point z on the circle from the 
inside (respectively outside). 

Next, observe the following factorization of the jump condition (|75|1 : 



(76) 



N!^(e'^) = N'!(e''') 








^{0) 








To take advantage of this factorization, we introduce two new contours T,± which together with E bound 
two concentric annular domains A± as shown in Figure |H1 




/\S_ 



Figure 8. The annular domains A± of the equivalent d problem for polynomials orthogonal 
on the unit circle. For a given e > 0, the contour T.+ corresponds to \z\ = 2^*^ and the 
contour I]_ corresponds to \z\ ~¥ . 
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We will now need some extension of the function (f>{d)^^ defined for z on the unit circle with arg(z) = 9 
to the annular domain A^ U A-. To make use of the family of extensions defined in |SJ), we now make the 
following assumption about the weight 0(0). 

Assumption 3.2. The function V is of class C''~^{S^) for some k ~ 1,2,3, ... . 

Note that when k = 1 this assumption is contained in Assumption 13. 1 1 but when fc > 1 it provides new 
information. Recall the "bump" function B with the properties listed in § 11.51 Then, for any integer m in 
the range 1 < tti < fc, and for any e > 0, we may apply the extension operator Em to the function V and 
therefore define a matrix P^ ei^^ ^) ^^ follows. 



(77) 



eir,0) 



1 







N"(z) 



N"(z) 



N"(z), 



z-"B(log(r)/e)e^'"^('-^^) 1 



-z"B(log(r)/e)e-E'"^(''-«) 1 



for z = re'^ G A- 



for z ~ re'^ G A 



+ > 



for z 



A9 



(^ A+\JA- 



^ from the unit 



Thus, the factor e-^''»^('''^) appearing above is our selected extension of the fmiction 
circle to the regions A+ and A_. 

Unlike M"(z) and hence N"(z), the matrix PJ^ ^(r, 6) is not piecewise analytic because the factors relating 
PJJj f_{r,0) to N"(z) in the domains A± are not analytic. Indeed, in view of lO, the exponent i*^,„\^(r, 6*) is 
not an analytic function. Note however that it follows from Assumption 13.21 and the analyticity of M"(z) 
and hence of N"(2:) for \z\ ^ 1 that the matrix PJ^ ^{r, 0) is continuous for ?- ^ 1 as long as 1 < tti < fc. In 
particular, the "bump" function factor i3(log(r)/e) ensures that P5^g(r, 6') is continuous across the circles 
S±. At the circle of discontinuity, E, the boundary values taken by P"j ^(r, 9) satisfy the jump condition 







(78) 



limP;j,,,(r,0) 

rjl 










<i Remark: Note that the approach to this problem taken in [H], where analyticity of (f) is assumed, 
amounts to replacing EmV{r,9) with the analytic extension EaoV{r,9) and omitting the "bump" function 
factor i?(log(r)/e), with the latter being at the cost of an exponentially near-identity jump discontinuity 
across the inner and outer circles S±. \> 

Next, we may remove the jump discontinuity along the unit circle by introducing a model matrix P(2;) that 
is analytic for \z\ ^ 1, tends to the identity matrix I as z — > cx), and that takes on the unit circle continuous 
boundary values P+(2) (respectively P_(z)) from the inside (respectively outside) that are related by 



(79) 



P+(e'^) = P_(e'^) 







Such a matrix can be found in closed form. Setting 



(80) 



Q(Z)=P(2) 



-1 



1 



for l^l < 1 and Q(z) = P(z) for \z\ > 1, one may equivalently seek a matrix Q(z) that is analytic for \z\ ^ 1, 
tends to the identity matrix as z — > oo, and that takes on the unit circle continuous boundary values Q_|-(z) 
(respectively Q_(z)) from the inside (respectively outside) that are related by the diagonal jump condition 



(81) 



Q+(e"') 



Q_(e^^)0(0)^^ 
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where a^ denotes the Pauh matrix 
(82) 



^3 




Clearly, we may seek Q{z) as a diagonal matrix. Assumption 13.11 guarantees that log{(f){9)) = —V{9) is 
a well-defined on the circle that satisfies a uniform Holder continuity condition with exponent v G (0; l]i 
and therefore we may obtain a matrix Q(z) with the aforementioned properties in the explicit form Q(z) = 
S^{zY^ where S^{z) is the Szego function associated with the weight 0(0) as defined in ||2SJ). Going back 
to P(z), the model matrix we will use to remove the jump discontinuity in PJJj ^{r, 9) for r = 1 is defined by 
the explicit formula: 



(83) 



P(z) := <^ 




UI>1 



\z\<l 



To actually remove the discontinuity, we introduce a new matrix function HJJj ^(r, 9) defined for r y^ 1 by 
the formula 



(84) 



H" (r,0):-P;:,Jr,0)P(z)-i. 



By Assumption 13.11 and Assumption 13.21 the matrix H"^ ^ (r, 6) is continuous throughout the two regions 
r < 1 and r > 1. Moreover, a continuous extension to r = 1 is possible because PJJ^ <:(?', 6*) and P{z) satisfy 
the same jump condition at r = \z\ = 1. Thus, we see that for 1 < m < k, and for any e > 0, the matrix 
function HJJ^ ^ (r, 9) defined by H84|l may be viewed as a continuous function on the whole plane with polar 
coordinates — tt < 9 < tt and < r < €«. 

At this point; we can summarize the explicit transformations we have introduced and relate H^ ^(r,6) 
directly back to M"(z). Combining ^^, {77)), ^^, and ^^, we have by definition 



(85) 



'aiJr,9) 



M"(z) 



M"(z) 



-S4,iz) 

S^z)-^ z"50(z)B(log(r)/e)e^'"^('-^'') 

z-"5^(2)-i 

54z)-iB(log(r)/e)e^™^('^'«) z^S^z) 



<r < 1 



r > 1 



where z — re'^. 

For r < 2~'^ and r > 2*^ we have i3(log(r)/e) = 0, and in these regions the matrix HJ^ ^.(r, 9) clearly inherits 
analyticity from M"(z); in other words, in these regions HJJj ^(r, 9) is a smooth function of the combination 
z = re*^. However, for z e Q,±, the matrix 'iii'^^{r,9) is certainly not analytic. In order to measure the 
deviation from analyticity in the case when m = k, we introduce a further assumption on 

Assumption 3.3. 

tinuous. 



The function V is of class C^ ^'^{S^), that is, the function V^^ {S) is Lipschitz con- 



note that since fc > 1, this condition implies in particular that V{9) is Lipschitz, and thus the Holder 
continuity part of Assumption 13.11 is subsumed. With Assumption 13.31 in force, we may compute the d- 
dcrivative of HJJ^ ^(r, 9) for r ^1 and for all integer m in the range 1 < m < fc. Differentiation of 18 5|) yields 
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(86) 



dU^4^,0) 



M"(z) 



M"(z) 



< r < 1 



r > 1, 



^7,^Jr,^ 





z"5^(z)9 [B(log(r)/e)e^-^('-'«)] 


for almost all 6, and then elimination of M"(z) in terms of if^ ^{r, 9) using H85|l again gives 
(87) dU'^Ar, 0) = H;^,,,(r, 0)W'^Jr, 9) , for r ^ 1 and almost all 9 e S^ , 

where the matrix W"^ ^ (r, 6) is given for r ^ 1 and almost all 6 by the explicit formula 

z"5^(z)2a [S(log(r)/e)e^'"^('-'^)] 





z-^'S^zy^d [B(log(r)/e)e^'"^('-''')] 
In particular, we see that the matrix HJ^ ^(r, 9) is a solution of the following 9-problem. 

d Problem 3.1. Find a 2 x 2 matrix \J{r,9) with the properties: 

Smoothness. XJ{r,9) is a Lipschitz continuous function throughout M.'^ . 
Deviation Prom Analyticity. The relation 

(89) dlJir,9)=\Jir,9)-W;^Jr,9) 

holds for all points in M" with the exception of a set of measure zero. The matrix WJJj ^ (r, 9) is 
defined almost everywhere by (|88|l and is essentially compactly supported. 
Normalization. The matrix U(r, 9) is normalized at r = oo as follows: 

(90) lim U(r, 9)^1. 



:= < 



< r < 1 



r > 1 . 



In writing down this 9-problem, we have focused on just a few specific properties of the matrix HJ^ ^(r, 9). 
However, it is important that in doing so, we have not introduced any spurious solutions. 

Proposition 3.1. Suppose that (p = e~^ where V : S^ ^ S. is of class C''^^'^{S^) for some fc = 1, 2, 3, ... . 
Then for all n = 0, 1, 2, 3, . . . , for m = 1, 2, . . . , fc, and for all e > 0, the matrix W^ ^(r, 9) is well-defined 
almost everywhere by l|88() and d Problem \S.l\ has a unique solution, namely U(r, 0) = HJ^^(r, 0). 

Proof. The existence of a solution follows from H85|) and the existence of M"(z) for n = 0, 1, 2, ... . To estab- 
lish the uniqueness we first consider the determinant of any solution of d Problem l3.1l Clearly, det(U(r, 9)) 
is a Lipschitz continuous function that tends to 1 as z -^ cx). Moreover, the relation 9det(U(r, 6*)) = 
tr(W;^,,(r,6'))det(U(r-,6')) holds almost everywhere, and thus by ^ we see that adet(U(r,6')) = holds 
almost everywhere in the plane. It follows that det(U(r, 6*)) is not only Lipschitz continuous, but is in fact 
an entire function of z = re that tends to 1 as z — > cxd. Therefore from Liouville's Theorem we see that 
det(U(r, 9)) = 1. Next, consider the matrix ratio of any two solutions U(7', 9) and U(r, 9) of d Problem 13. II 
this is the matrix R(r, 9) defined by 



(91) 



R(r, 



V{r,9)V{r,9y 
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Since det{\J{r,9)) = 1, it follows that R(r, 0) is Lipschitz continuous throughout the plane. By direct 
calculation, we have 

(92) = U(r,0)W;;,_,(r,0)U(r,e)-i-U(r,0)W;j,,,(r,0)U(r,0)-i 

= 0, 

holding almost everywhere in the plane. It follows that R(r, 9) is an entire function of z = re'^ that tends 
to the identity matrix as z — » cxd, so again by Liouville's Theorem we get R(r, 9) = I, or equivalently 
tS{r,9)=V{r,9). D 

The unique solution of d Problem 13.11 can also be expressed as a solution of an integral equation with 
Cauchy kernel. 

Proposition 3.2. Suppose that cfi = e~^ where V : S^ ^ M. is of class C'^~^'^(S'^) for some k = 1, 2, 3, ... . 
Then for all n = 0, 1, 2, 3, ... , for m = 1,2, ... ,k, and for all e > 0, the matrix W^ ^(r, 9) is well-defined 
almost everywhere by (|88|l and the corresponding solution XJ(r,9) ~ HJ^ e(''i^) '^f '^ Problem \8.1\ satisfies the 
integral equation 

TT J J Z' - Z 

where z = re^^ , z' = r'e*^ , and dA' is a positive area element dA' — r' dr' d9' . The integral is taken over 
the entire plane. 

Proof. Recall that the Cauchy kernel is a fundamental solution for the d operator. In the relation H87|l we 
may replace dti'^^{r,9) by d\H.'^^{r,9) — I]; multiplying by the Cauchy kernel and integrating over the 
whole plane gives the identity 

(94) 1 I I 9[H;^,,(r', 9') - I] ^^^, - ^ f f ^".^(^'' ^'W'^nAr', 9') ^^^, 

n J J z' - Z TT J J z' - z 

On the left-hand side the d operator differentiates with respect to the primed variables. Since the Cauchy 
kernel is absolutely intcgrable, we may evaluate the integral on the left-hand side by replacing the domain 
of integration by the region \z' — z\ > 6 > and subsequently taking the limit S ^ 0. For each positive S we 
may apply Stokes' Theorem and use the facts that d[{z' — z)~^] = for \z' — z\ > S and that H"^ ^{r, 9) —I 
tends to zero as r — > oo to evaluate the integral over the region \z' — z\ > (5 in terms of a line integral over 
the boundary. Thus wc have 

sio 2tt5 J|2'_^|=i- ' TT J J z' ~ z 

where d(,' is an arc length element. From the continuity of H"jj(r, 6') the integral equation (|93|l with 
U(r, 9) = H;^^,(r, 9) follows. ' D 

3.2.2. Asymptotic solution of the integral equation. Estimates of H.^ ^{r,9) and its derivatives for large n. 
Being as knowledge of the matrix H"j ^{r, 9) is equivalent to knowledge of M"(z) and hence of the polynomial 
of degree n in the system of polynomials orthogonal on the circle with respect to 0, wc would like to use the 
integral equation H93|) to characterize H"j ^(r, 9). There is a difficulty in that, while existence of solutions for 
()93|l is not an issue, one does not automatically have uniqueness. However, it turns out that if the parameter 
n is sufficiently large, then the integral equation H93|) defines a contraction mapping and thus may be solved 
by iteration yielding a unique solution in the form of a Neumann series. In this connection, we can also 
obtain from (|93|l asymptotic information about the matrix HJJj ^{r,9), and consequently of the orthogonal 
polynomial 7r„(z), in the limit n -^ cx). 

In order to study H93|) , it is useful to characterize the family of matrix functions WJ^ ^ (r, 9) more concretely. 
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Proposition 3.3. Suppose that V : S^ ^ 'R is a real function of class C'^^^'^{S^) for some k > 1, that 
m is an integer satisfying 1 < m < k, and that e > is fixed. Let the integer D he defined as D := 
niin(fc — m,m — 1). Then, the matrix function W"^ ^ is of class Cq ~ ' (M^ \ {0}) if D > 0, and of class 
I/g°(R^ \ {0}) if D ^0. Moreover, if a and (3 are nonnegative integers such that a + f] < D, then there is a 
constant Cmje > such that for all n the estimate 



(96) 



Qa+P 



■W:j,ir,0) 



<C(,"f)n''e-"l'°s('-)l|log(r)r 



^nP|log(r)|P 

p=0 



dr'^dOP 
holds throughout the region |log(r)| < elog(2) containing the essential support ofW"-^^(r,9). 

Proof. WJJj g(r, 9) vanishes identically outside of the annulus | log(?')| < elog(2). In the disjoint regions r < 1 
and r > 1, the matrix function WJJ^ ^(r, 6) is infinitely differentiable with respect to r, and the issue is the 
continuity of these derivatives at r = 1. The relation {TJ implies that for each fixed 6, W"^ ^(r, 6) is propor- 
tional to (log(r))™~-'^ near r = 1 (where B(log(r)/e) = 1 holds), and thus all derivatives of WJ,\ ^(r, 9) with 
respect to r through order m — 2 are Lipschitz continuous at r = 1, and the derivative d"^~^W^ ^{r, 0)/dr"^~^ 
remains bounded as r ^ 1, but experiences a jump discontinuity at r = 1. 

On the other hand, if r ^ 1 is fixed, then from (jSJ and iQ, the matrix W"j^(r, 6*) depends analytically 
on derivatives V'^^^O) for < j < to. Since V is of class C''^^'^(S'^), all derivatives of WJ^J^ ^ with respect to 
9 through order k — l — m will be Lipschitz continuous, while the derivative 9'^~™WJJj ^(r, 9)/d9^~'^ will be 
defined for almost all 9 and will be uniformly bounded. 

To have all mixed partial derivatives of total order at most D — 1 to be Lipschitz continuous, it is sufficient 
to have both D <m~l and D < k — m. If for 1 < to < A: these ineciualities force D = 0, then no derivatives 
of WJJj J (r, 6') may be taken at all, but WJ^g(r, 0) is uniformly bounded and compactly supported in the 



annulus | log(r)| < elog(2), that is, W^J, , G L^{W \ {0}). li D = min(fc 



1) > 0, then we learn 



that W" P C, 



D~l,l, 



\{o}). 



Now there are absolute constants C± > such that 



(97) 



ga+0 



dr^d9P 



W" ir 



C4 



^ W„,6,a",/3"(r,( 



OL'+(x"=Oi 

0'+l3"=l3 



Qa'+P' 



jz'^S^izf 



r < 1 



C- 



J2 ^" 



.fi"{r, 



^ dr°''d9P 



jz-''s^{zy 



P'+f3"=l3 



r > 1 



where 
(98) 



Wm,e,a",P"{r,9) : = 



aa"+P" _ 
d 



S(log(r)/e)e^™^('^'^) 



dr°"'d9f^" 

Generally, the derivatives indexed by a' and /?' are uniformly bounded throughout the regions 2~'^ < r < 1 
and 1 < r < 2*^ by a constant multiple of n" ^'^ e~"l '"^^''■'1. Furthermore, using JTI), the derivative 
Wm,e,a",i3"{r,9) is Uniformly bounded by a constant multiple of | log(r)|"'^^^" throughout the region 
I log(r)| < elog(2). Therefore setting a" = a — a' , an inequality of the form 



(99) 



ga+p 



dr"d90 



w:;,,(r,( 



<C'i"f)e-"l'°sMI J2 ^7i'3'|log(r)r"i~"n"'|log(r)|"' 

a'=0/3'=0 

holds in the region | log(r)| < elog(2), where Cm,e > is a constant. 

Finally, since nf^' < n^ , the inequality ^ follows, where C^:P = (/3 + l)^^;/^ 



D 



An important part of our analysis will be the estimation of certain two-dimensional Laplace-type integrals 
with Cauchy kernels. The main workhorse in this connection is the following lemma. 
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Lemma 3.1. Let e > and v >1 he fixed constants. Then there exists a corresponding constant K^^, > 
such that the estimate (note that z' — r'e*^ ) 



(100) 



sup 



2' 



r'dr'e-"l'°«('''"|log(r')r"^ 



d0' 



<K, 



log(n) 



holds for sufficiently large n. Moreover, for each p > 2'^ , there exists a constant K^,^p > such that the 
estimate 

r2' 



de' 



K, 



€.,l',P 



(101) sup / r'dr'e-"l'°s('-)l|log(r')r"M 7-> T^ 

holds for sufficiently large n. 

Proof. As 9' varies over S^, the minimum value of |z' — z| is achieved at 0' — 9, and we thus have \z' — z\ > 
\r' — r\, which imphcs the inequahty 

(102) r^^< 2^^ 



Let // be a positive constant. Clearly, there is another positive constant Ci depending on /i but not on r or 
r' such that 



(103) 



n r - ^1 

^TT W ~ Z\ 



because the condition | log(r') — log(r)| > fj. also bounds \r' — r\ away from zero. Furthermore, if fi is 

sufficiently small, there is a positive constant C2 depending on fi but not on r or r' such that 

(104) 

— < C2 log n + - — — — - 

_^\z'~z\ \ |log(r')-log(r)|. 

To establish (|104|l . note that for /i sufficiently small the condition | log(r') — log(r)| < /i implies that \r' — r\ < 
tt/2. Assuming without loss of generality that 6* = 0, we use the following estimates of the integrand. For 
\9'\ < \r' —r\ we use the estimate \z' — z\ > \r' — r\ which follows from the triangle inequality applied to the 
identity z + {z' - z) = z'. For \r' - r\ < \9'\ < tt/2 we use the estimate \z' - z\ > r'\ sin(6'')| > r\9'\/2 which 
follows from the diagram shown in Figure|^ Finally, for tt/2 < \9'\ < tt we use the estimate \z' — z\ > r which 





1 


= r'e''' 






r' / 




r' sin(^') 


z' ~ 


Z/ 


/e' 




n 




^"^\^ 



z = r 



Figure 9. The estimate \z' ~z\ > r'\ sin(0')|. The right-hand side can be replaced by r'\9'\/2 
for \9'\ < tt/2. 

follows from the Law of Cosines because cos(6'') < 0: \z' — zp = r^ + (r')^ — 2rr' cos{9') > r^ + (r')^ > r^ 
(again, see Figure EJ- Combining these estimates, we have 



(105) 



d9' 



4 r/2 d9' TT 
<2 + - ^^ + - < 2 

\z' - z\ - r' ./i^,_^i 9' r - 



22+^ / Z- + 2^7re^ . 



This is clearly bounded above by a constant multiple of logd?'' — r\^^) for /i sufficiently small, and from 
this (|104|l follows as well (adding a constant inside the logarithm keeps the bound positive away from the 
singularity as long as n > 1, and for later purposes it is convenient to take the additive constant to be n). 
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Now we estimate the integral over r' . Using r' < 2*^, along with H1U3|) and (|104|) . and changing the 
integration variable from r' to s = 7ilog(r'), we get 

'' r'dr'e-"l'°«('-')l|log(r')r-^ T ^^ < ^-^ [ dse-^'^lsr^ 

J^^ \z' - z\ n-" Jj^ 



(106) 

where 
(107) 

Finally, we have 
(108) 



2^^C. 



Jb 



nlog(r)| 



A := {s such that \s\ < rielog(2) and \s — nlog(r)| > nfi} 
B := {s such that \s\ < 7ielog(2) and \s — nlog(r)| < nij,} . 

/•oo 
A 



which is finite and independent of n and r because v > 1, while 

n 



I dse-\%\''-^ (n + 



s — rtlog(r)| 



< 



ds e 



log n 



|s-nlog(r)| 



(109) 



log(n) 



dse-l^l|sr-i 



dse-l^llsr-i 



log 1 



1 



and the last integral is bounded independently of n and r, since by Cauchy-Schwarz, 
(110) 



dse-l-'l|sr-Mog(l 



nlog(r)| 



< 



dse-2|-^IU|2-2 



1/2 



nlog(r)| 



1 



ds log 1 + 77 



-,1/2 



with both factors being finite. Thus, an upper bound for the integral of interest is proportional to log(n)/n'' 
in general, which proves pOOfl . If p > 2'^ and |log(|z|)| > log(p), then it is not necessary to divide the 
integration into sets A and B, and the bound Ij 103(1 can be used over the whole range of integration in which 
case the upper bound is then proportional to \/ny which proves ((101(1 . D 

With these results in hand, we can formulate and prove the following. 

Proposition 3.4. Suppose that (p = e^^ where V : S^ ^ M. is of class C''^^'^{S^) for some k ~ 1, 2, 3, ... . 
Let the integer m lie in the range 1 < m < k and fix e > 0. Define the integer D := min(/s — to, tti — 1) > 0. 
Then, for all n > the matrix WJ^ ^(r, 6) is well-defined almost everywhere by 1(88(1 . and for all n sufficiently 
large, HJJj ^(1^,9) is given by a Neumann series 

(111) h;^,.(^, ^) = I + (>vr„,ei)(r, e) + (w;;,, o >v,"„,,i)(r, 0) + --- 

which converges in the norm \\\ ■ \\\d, where the double integral operator W," ^ is defined by 

F{r',e')W^Jr',P'^ 



(112) 



(W,"„,,F)(r,0):=-i||. 



■dA' . 



In particular, if D = then HJ^ ^ lies in the space L°°(M.'^), and if D > then HJ^ ^ lies in the space 
C^~^'^(M.'^) and \\\tiln cIII-D "is finite. For all integer p in the range < p < D, the following estimates hold 
for sufficiently large n: 



(113) 



IH^ 



< ciS\ 



logjn) 
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(114) 



|H" -I-W" I|||p<CiP) 



^ j-,2i7i — 2p 



^(p) 



where Cm'e > is a constant. Furthermore, for each p > 2*^ and for all integer p in the range < p < D, 
the following estimates hold for sufficiently large n: 



(115) 



(116) 



E 

a+l3<p 



sup 

-7r<e<7r 
log(r)|>log(p) 



dx'^dyP 



[H",(r,( 



< ci^l 



m,p ^m-p ' 



E 

a+l3<p 



sup 

-7r<e<7r 
log('")l>log(/o) 



Qa+0 



dx'^dy'^ 



[H;^,,(r,0)-I-(W,"„,eI)(r,0)] 



< ^(P)2 log(") 



™>P ^2m-2p 



=.(P) 



where Cm.p > is a constant. 

Proof. Fix fc > 1 and rn in the range 1 < m < k, and set D = min(fc ~ m,m ~ 1). Let p be a nonnegative 
integer satisfying p < D. If p > 0, suppose that F{r,6) is matrix function of class Cp~^'^(M^) with all 
derivatives of total order no greater than p uniformly bounded in the whole plane, and if p = 0, suppose that 
F{r,6) is of class L°°{M.'^). For such F(r,6) we recall the norm (I15II . where the cartesian coordinates x and 
y are connected to the polar coordinates r and 9 in the usual way: x — rcos{9) and y ~ rsm{6). Because 
for z' ^ z, 



(117) 



aa+0 1 fia+0 



dx°'dy'^ 



dx'^'dy'P 



for a > and /3 > with a + /3 < p, we have 



(118) 



dx°'dyP 



(W,"„,,F)(r,0) 



(_l)a+/3 



F(r',0')W,"„,,(r',0' 






dx'-^dy'P 



dA' 



f)a+0 JAI 

[F(r',0')Wl,(r',0')]-r 



TT J J dx'^dy'^ " ■ ■ """ ' " z' — z 

Note that in order to integrate by parts in 1118|l for all a and /? of interest, we must have p < k — m. Now, 
since WJJ^ ^(r, 0) is compactly supported in the annulus | log(r)| < elog(2), 



Qa+fi 



dx'^dyf^ 



[F(r,0)W:;,,(r,0)] 



< 



E 



a.' -\-a"^Q. 
/3'+f3"=l3 



sup 



ga'+/3' 



Fir,d) 



(119) 



< iiiFiiip-E E 

q"=0/3"=0 



dx'^'dyf^ 

Qa"+P" 



Qa"+P" 



dx""dy' 



77^1. Ar,d) 



dx°'"dyP 



-W" (r 



7n,e V ' 



< iiiFiii, Y. ^r:'/ 

a'+l3'<a+l3 



Qa'+P' 



QrO^'dOP 



T^lAr.O) 



'("',/3') 



where Km,e are some positive constants. Using Proposition l3.3l we then find that 



(120) 



Qa+I3 



dx'^dyP 



[F(r,0)W:;,,(r,( 



Q+/3 



< Kl^y^^ ■ |||F|||pe-"l'°s('-)l|log(r)r"-i-("+'^^) ^ n-''|log(r)P 



,(0 



where Km,e > is some constant. Finally, we arrive at the estimate 

f>{p) P /•2' 

(121) |||>v;;,,F||| < " 
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(-IPJ P 1-2- i-TT ini 

^^^•|||F|llp-E"' '/ ^''^^'e""l'°^^"'^l|log(r')r+^'-^-M ^P- 



Hp) 

■m,e 

Set v := m + i — p, and note that v >1 since as j ranges from Q to p, v ranges from m — p to m, and we 
have p < D < m — 1. Lemma mi may thus be apphed to each integral on the right-hand side of H121|l . with 
the resuh that 

(122) |||W™.F||| <^2^ii^iM|||F|||,, 

for some constant Cm,e > and n sufficiently large. We note in passing that in order for (|122|l to provide 
control of the operator W^ j, we need to have p < m. The two restrictions in force on p, namely p < m and 
p <k — m, have been expressed in the statement of the proposition as the inequality p < D. If one restricts 
attention to those z = re'^ for which | log(r)| > log(p) > elog(2), then Lemma [3.11 implies the inequality 



(123) Y. ^^P 



|logWI>log(p) 






dx'^dyP 



< 



^m—p ' 



where Cm'p > is a constant and n is sufficiently large. 

From (|122|l it is clear that if n is sufficiently large, the double integral operator W^ ^ defined by the formula 
(|112|1 and acting in the integral equation (|93|1 thus defines a contraction mapping in the space C'~^'^{M.'^) 
equipped with the norm ||| • |||p, or in the space L°°(R^) Up = 0. This implies that there is a unique solution 
of H93|) in this space that may be found by iteration resulting in the ||| • |||p-convergent Neumann series 

(124) U(r, e):^l+ (W;; J)(r, 0) + (W;;,, o W:^J)ir, 0) + ■ ■ ■ . 

In particular, choosing p — one sees that the Neumann series (|124|l furnishes a unique solution of the 
integral equation (|93|l in the space L°°{M.'^). Since Proposition 13. 21 guarantees that the matrix HJ^ e('':^) is 
a known solution of the integral equation (|93|l that is (in particular) uniformly bounded in the plane, we 
may identify it with the Neumann series H124(l for n sufficiently large. From this point forward in our proof, 
we assume that n is indeed large enough for this to be the case. Since the same Neumann series H124|l also 
converges in the norm ||| • |||p where p may be taken to be as large as D, we also learn that if 13 > then 
Hl,^,{r,e) lies in the space C^-^'\R'^) and that |!|H;^„,J||c is finite. 
From ()124|l . taking ||| • |||p norms and using (|122|l . we see that 



(125) mi,-mp<mT.( nZ ] ^^ 



C ff.l0g(n) y ^ Cff.log(n) 



if n is large enough that n™ ^ exceeds 2Cmi log(n). Since according to the integral equation (PS|l satisfied 
by H.'^ ^(r, 0) we have 

(126) ll^^,ir,0)-l-{W:i,J){r,0) = (W,"„,,(H;^,, - I))(r, 0) 

we may take norms and use (|122|l and H125|) to learn that 

(127) m-^,, - 1 - w:,^Mp < m\ ";;L-2; 

holds for sufficiently large n. The proof of the estimates (|113|l and (|114(l is complete upon appropriate 

(p) 
redefinition of the constant Cm'e- 

Note that p23|l implies that the upper bounds in (|125() and (|127ll can be reduced by a factor of log(n) if in 

each case the supremum on the left-hand side is taken over only those values of r satisfying | log(r)| > log(p) 

for a fixed p > 2*^. This completes the proof of the estimates (|115|l and (|116|l upon appropriate redefinition 

of the constant Cm'p- □ 

< Remark: To uniformly control p derivatives of H"^ ^(r, 9), Proposition l3.4l reauires that m should lie in 
the range 1 +p < m < k—p, and therefore in order for there to exist suitable values oi m, V : S^ -^ M. should 
be of class C''~^''^{S^) for some fc > 2p + 1. Also, note that the utility of the estimates H114|) and Hlltj|l is 
that the matrix WJJ^ ^(r, 9) is off-diagonal, so the diagonal matrix elements of H'^ ^{r, 9) — I experience more 
rapid decay than do the off-diagonal elements. > 
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3.3. Proofs of theorems stated in § 13.11 If we solve (|85|l for M"(z) in terms of HJJj e(?', 6*), then since 
Mii{z) = 7r„(z), the monic polynomial of degree n in the system of polynomials orthogonal with respect to 
the inner product (•,•)</, defined by 0, we can easily obtain from Proposition 13.41 asymptotic formulae for 
7r„(z) and its derivatives, valid for large n, with uniform error estimates. 

3.3.1. Asymptotic behavior o/7r„(z) for \z\ > 1. Proof of Theorem \3.1\ In the region r > 1 we have for each 
e > and for each m = 1 , . . . , fc the exact representation 

(128) ^„(z)=if;;,,,^i,(r,0)z"e^(^)-7j;;,^,^i2(r,0)i?(log(r)/e)e^'"^('-''^)-^(^). 
Here we have used H38|l to write S^{z) = e^*-^^ for \z\ > 1. Equivalently, 

(129) ^„(z)z-"e-^(^) - 1 = [i7rn,Mi(^' ^) - 1] " ^-"^™,M2(^ 0)B(log(r)/e)e^'"^('-'«)-2^(^) , r > 1 . 
If p > 1 is fixed, then we may choose e > small enough that i3(log(r)/e) = whenever \z\ > p. In this 
case, we have simply 

(130) ^„(z)z-"e-^W-l = if;;,^,^ii(r,0)-l, \z\ > p > 1 . 

The best decay estimate comes from taking m — k. In this case, using Proposition |^31 (specifically recalling 
the estimate (|116ll and the fact that WJ^ c(^:^) is an off-diagonal matrix), we see that for some constant 
ifp > 0, 



(131) 



sup |i7" .ii(^^)-l|<i^, 



-■7r<e<ir 
r>p>l 






Now as the combination 7r„(z)z "e ^'^^ — 1 is a function of z that is analytic in the region \z\ > 1 and 
decaying to zero as z -^ cx), we may express arbitrary derivatives of it as Cauchy integrals: 



(132) 



dzP 



7r„(z)z-"e-^(-) - 1 



^ni 



7r„(s)s" 



-N{s 



1 



■ ds . 



where the contour of integration is oriented in the counterclockwise direction, and \z\ > p. Using (|131|l to 
bound the integrand then gives a uniform bound of the same order of magnitude for derivatives over regions 
bounded away from the circle \z\ ~ p, which can be taken arbitrarily close to the unit circle. This proves 
(|40|) . and completes the proof of Theorem 13. II 

3.3.2. Asymptotic behavior of iTniz) for \z\ = 1. Proof of Theorem \3.'A Theorem 13.21 follows from the 
estimate H42|l by noting that the error is an analytic function of z in the exterior domain |z| > 1 that decays 
as z — s- oo, and therefore ()42|l implies the more general result stated in Theorem 13.21 via the maximum 
modulus principle. To prove (|42|l . we fix any positive value of e and consider 1 < r < 2*^'^ in which case 
i3(log(r)/e) = 1 and therefore (jl29|l implies the following formula: 
(133) 
dP 



dzP 



7r„(z)z-"e-^('~^ ^l =dP K,,,n(^e) - l] -d" z-H:^^,^,^{r,9)e''-''^^''^-^''^^^ 



1 < r < 2 



,/2 



We remind the reader that estimates on derivatives like p33f) are valid forl+p<m<fc— p (see the 
remark at the end of i) l3.2() . Using the estimate (|114ll from Proposition 13 . 41 and noting that WJJj ^(r, 6) is an 
off-diagonal matrix, we see that for some constant Kp > 0, 

\og{nf 



(134) 



snv\dP[H:^^,^^,{r,e)-l]\<K, 



P„2m-2p ' 



if r7i < k ^ p. On the other hand, the dominant contributions actually come from those terms in the second 
member of the right-hand side of (|133|l in which none of the p derivatives fall on the exponential factor 
^E^v{r,9)-2N{z) (^■^}^[(.\y j^^s k — m+ 1 uniformly bounded derivatives). Since \z\ > 1, it suffices to estimate 
n^dP~^ H"^ J 12 (''i^) with the use of the inequality l|113() in Proposition l3.4l Taking m = fc — p for the best 
possible decay estimate then gives 



(135) 



sup 

-TT<e<TT 

l<r<2'/2 



dP 



dzP 



■^„(z)z-"e-^(^) - 1 



<K. 



log(») 

P ,-,fc-2p 
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where Kp > is a constant. This proves (|42|l . upon taking the Hmit r | 1 and writing the z derivatives in 
terms of 6 (differentiation commutes with the Hmit process). 

3.3.3. Asymptotic behavior o/ 7r„(z) for \z\ < 1 and of ^n ■ Proof of Theorem 15'. 51 Theorem \cl.4\ and 
Theorem \3.5l Using H85() and the fact (sec H38() ') that 5^(0) ~ e"^", we have the exact representation: 

(136) 7,1-1 = -M',\{0) = i7;;.,,,22(0, e)e^° , 

and, whenever \z\ < 1, 



(137) 7r„(z) = M-,{z) = z"e-^«-^(i/^)i?(log(r)/6)e^'"^('-'«)H,"n,.ai(^ 0) ™ e''°+''^'^^'>H:i,^,^,,{r, 9) , 

where (|38ll has been used, and z = re*^. 

To prove Theorem 13.41 we simply apply pit)|) from Proposition 13.41 in the case p = and m — k to the 
identity (|136|l . This immediately yields H5U|) and completes the proof. 

The proof of Theorem 13.31 is based on a similar analysis of p37f) . Recalling that p G (0,1), one can 
choose e > -21og(p)/log(2) and then i?(log(r)/e) = 1 in p37|) for p < \z\ < 1. The estimate (gHl) then 
follows by taking m = k, and using H114|l from Proposition 13.41 in the case p — 0. Similarly, choosing 
e < — log(p)/log(2), we have i3(log(r)/e) = in H137() for \z\ < p. Again taking m = k, one obtains H47() by 
using H116|l from Proposition 13.41 with p = 0. This completes the proof of Theorem 13. 31 

The rest of this section will be devoted to the proof of Theorem 13.51 Wc begin with H137|) for m ~ fc, 
a formula that is valid for all z with \z\ < 1. Using (|114|1 from Proposition 13.41 and keeping the term 
corresponding to VV^^I, we arrive at the formula 



n J J z' ~ z \ n^*^ / 



2-^<r'<l 



as n — > oo, where the error term is uniform for \z\ < 1. Here, dA' ~ r' dr' dO' is an area element, and in the 
integral 9' varies over S^ (the support of WJ! ^ ]^2(''' ^) i^ the annulus given by the inequalities 2^'^ <r<l). 
In the annulus of support of M^^^ ]^2(^' ^)' ^^ have from H88|) that 



(139) Wl,^^^{r,9) = d z"e-2^°-2iV(i/")B(log(r)/e)e^'=^(''^^) 

< Remark: The fact that WJ} ^ j^jl'^i^) i^ in the range of d means that the double integral in (|138|) can 
be reduced without approximation to the sum of an explicit contribution and a contour integral. Indeed, by 
the inversion of the d operator, 

1 /■ /■ M/" (r' 0'\ 

(140) -^ / / y_V '^^' = ^^^^ + z"e-^^''-^^(i/^)i3(log(r)/6)e^''^('--^)x(2-M)(0 , 

2-'=<r'<l 

where xi denotes the characteristic function of an interval /, and where F{z) is function analytic except on 
the circles \z\ = 2"*^ and \z\ = 1 bounding the support of WJ}^ i2('';^) that is chosen to make the right-hand 
side continuous and decaying as z — > cxd. These latter properties uniquely identify F{z) with the Cauchy 
integral 

1 f gU -Va+iO.(a,T^(s)) 

(141) F{z)^--—(i> ds, 

where the contour of integration is oriented counterclockwise. Note that it is the presence of the bump 
function _B(log(r)/e) that makes F{z) continuous at \z\ — 2""^. Unfortunately, this interesting formula, while 
apparently simpler than a double integral, is not as useful for asymptotic analysis as the alternative approach 
wc now follow. > 
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Continuing our analysis, we carry out the differentiation in (|139|l in the region |z| < 1 with the use of |(7J): 



(142) 



W"m,i2(^'^) = ^"e~2yo-2W(i/z)g£,y(r,e) ^B{log{r)/e) + Bilog{r)/e)dEkVir,e)] 



^n^-2Vo-2N(l/^) ^EkV(rfi) 



^e^'v(^\e), ^., _,„,,,_! 



aB(log(r)/6) + — ^^-^(-»log(r))'=-^i3(log(r)/e) 



and in the special case that r = \z\ > 2 ^/^, we have i?(log(r)/e) = 1 and dB{\og{r)/e) = 0, so 



V^^\9) 



(143) W^, ,,{r,e) = ^n^-^V„-2NilM^E,Vir^9) . 1^ ^_^ (_, log(^))fe-l ^ fo^ 2~'/' < r < 1. 

2r (fc — 1)! 
The presence of the z" factor together with the absolute integrability of the Cauchy kernel in two dimensions 
means that 

,144) -1 U ^jk^i^,A'^J- ff '''•^'■■-^'■•0\a' + O(2-"') 

TT J J Z' - Z TT J J Z' - Z 

2-'<r'<l 2-'/^<r'<l 

holds as n ^ oo uniformly for all z with \z\ < 1. Therefore using the simpler formula (|143|) in the integrand 
and integrating over the smaller annulus 2^"^ " < r' < 1 introduces an error that is uniformly exponentially 
small for \z\ < 1. 

Using the simple identity 

(145) d{- logir))'^ = -^{- logir))"-^ , 
we may rewrite (J143II in the form 

(146) Vt^,",.i2(r-,0) = -^z"e-2^"-2Ar(i/-)g£,y(ne)^(fe)(^)^(_i^g(^))fc ^ ^^^ 2-^/^ < r < 1 . 

Our subsequent analysis will be specialized to the case where V'^'^'^ (6) is piecewise continuous, with jump 
discontinuities at £ < oo angles — tt < ^i < • ■ • < 0^ < tt, and is (at first) only Lipschitz between the points 
of discontinuity. Then, between the points of discontinuity, y('^+i) (0) exists almost everywhere and may be 
identified with a bounded function. Under these circumstances, we may "integrate by parts" {i.e. apply 
Stokes' Theorem) with the following formula 

(147) 1 f f fir',e')dg{r\e') ^,^,_ 1 f f djfjr' ,e')gir' ,6')] ^,^, ^ 1 f f g{r' ,e')df{r\e') ^^^, ^ 

TT J J Z' - Z TT J J Z' - Z TT J J z' - Z 

and the first integral on the right-hand side may be exchanged for a sum of explicit terms and a contour 
integral as described in the above remark. In H147|l and in the rest of the proof, whenever the operator d 
appears in the integrand it acts on the primed variables. 

To prepare to use this technique, we begin with 
(148) 



1 [[ »w',^o,^,_.;f 

7' — ■/ ht ^-^ 



TT J J Z' — Z fc! ._ 

2-«/2<r'<l ^~^ 



1 /• /■ {z'Yh{r\&)V^^\Q')d{-\og(r')f ,., 

CLJ\ 



TT 

2-'/2<r'<l 



where /i(r, d) is shorthand for the following terms: 

(149) h{r, 9) := e-2^«-2lv(T7f)gS.v(ne) ^ 

and Ij refers to the interval in S^ of initial angle 9j and final angle equal to the point of next jump 
discontinuity as the circle is traversed in the counterclockwise direction. Using p47|) . we therefore find 

1 f f W"' ir' 6') 

(150) — / / ''''f^ ' ' dA' = -.h{r,e)+Kk{r,9), 

TT J J Z' - Z 

2~'/^<r'<l 
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where 



Mr,0) -.^ '-Y. 



i=i 



(151) 



1 

TT 



(zO"^[/^(r^gOFW(^0](-log(rO)fe ^^^, 



2-«/2<r'<l 

e'eij 



Kk(r,( 



fc! 



E 



1 

TT 



9 [(z')"/i(r', 6l')t^^'=H^')(- log(?-'))'' 



■dA' 



2-'/2<r'<l 

O'eij 



The double- mtegral expression Kkif, 0) may be reduced to contour integrals as follows: 

fc! 



(152) 
where 



Gfe(z) 



X,(r,0) = '-^z^h{r,e)V^''\e){~\og{r)fxi2-^n-.i){r) + Gu{z), 

s"Mkl,arg(s))yW(arg(s))(-log(|s|))^ 



i*^ 1 
fc! ' 27ri 



ds 



|s|=2-«/2 



S — Z 



(153) 



f aW 



+ fE 



fc! ^ 27ri J2-e/2^^Dj 

.7 = 1 



s"/i(|s|,0,)(-log(|s|))^ 



ds , 



s — z 



and A^ := V^'''^{9j+) — ^^'^-'(6'^— ). Since the ||| • |||o,i norm of the numerator in the first Cauchy integral 
is proportional to 2^^"/^, we may also write Gk{z) in the form 



(154) 



G'=W = iiE 






fc! ^ 2TTi Ja-^/^e^'j 



^MNI,^,)(-iog(|.|))^ 



s — z 



ds + 0{2 



-en/2\ 



as n — > oo where the exponentially small error term is uniform for \z\ < 1 (right up to the contour of 
integration; this is a consequence of the Plemelj-Privalov Theorem 122)). 
Let us now consider Jk [r, 9) . Note that 



d[h{r,9)V^''\e)] 



(155) 



V^''\9)dEkV{r,9) + ^V'^''+^\9) 



h{T,e) 



IZ 

2^ 



^(fc)(0)2 



(-ilog(r))^-i+y('=+i)(0) 



h{r,e). 



(fc-1)! 

Therefore, inserting this into the integrand for Jk{r, 9) and applying Lemma l3. II with i^ = 2fc to the integrals 

resulting from the first term above, we find 

(156) 



Jkir,9) 



Ak+l 



(fc + l)!E 
'log(n) 



1 

TT 



2-'/2<r'<l 



-^^(-log(-'))'= 



{z')''hir',9')V^''+^\9') 

■ (jyi 



O 



-,2k 



where the error is uniformly small as ri -^ cx) for all \z\ < 1. With y('^+^)(6') bounded in each Ij, we could 
in principle apply Lemma |3. II to the remaining integrals. However, this would only give a bound of order 
log(n)/n'^+^, and as Gk{z) will turn out to be (for most z) of size l/n*''+^, and we will want to consider 
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Gk{z) to provide the dominant term, we need to impose additional conditions on T/C^+i) in each Ij to see 
that Jk{r,9) is indeed subdominant. Therefore, we first use (|145|l to write Jk{r,6) in the form 



(157) Jkir,e) 



,-fc+i 



(fe + 1)!^ 

.7 = 1 



1 

TT 



(z')"/i(r', e')V^''+^'> {e')d{log{r'))''+^ 



■dA' 



2-»/2<r'<l 



O 



fog(n) 



,2fc 



The double integral above is of the same form as the original double integral we arc trying to compute (see 
(I148|l '). but with k replaced hy k + 1 (everwhere it appears explicitly; in the function h{r, 0), k remains k). 

To continue the analysis of Jk {r, 9) , we therefore use the further assumption that in each each of the 
intervals /_,-, V^(''+^'(6'') is a Lipschitz continuous function so that y('^+2)(6*') exists almost everywhere in Ij 
and can be identified there with a bounded function. Repeating the above steps, we find that 

(158) Mr, 6) = -Jk+i{r, 0) + A'fc+i(r, 0) , 

and with the help of the identity 

"yW(6))F(''+i) 



(159) 



d[h{r,0)V^'^+'\9)] = — 



-{-i\og{r)f-^+V^''+^\0) 



h{r, 9) 



(fc-1)! 

we may apply Lemma 13.11 to the integrals that result from substituting the above into Jk+i{r,6) with 
v = 2k + I (from the first term above) and with v = k + 2 (from the second term above). Consequently, 

'log(n) 



(160) 



Jk+iir,e)=0 



,fe+2 



where the error is uniformly small as n -^ oo for all r < 1. As before, /Vfc+i(r, 0) is given by (|152|l with 
Gk+i{z) given by (|154|l (note that in substituting fc + 1 for fc in these formulae, one leaves h{r, 6) alone). 

Combining these results, we have shown that 
(161) 



1 

TT 



wr,,{r',0') 



dA' 



'/2<r'<l 



k\ ^-^ 2m 

j=i 



2-'i/2e*''j 



(fe+1) 



(fc 



jfe+1 ' A^^ ' /■ 
VTV)\ j^^ 2m X 



5»Mkl,e.)(-iog(kl))'=+^ 



■ ds 



2-'/2e'"j 



o 



/log(; 



5l)+0(r"|log(r)|%-./.,i)(r)) 



holds uniformly for r = \z\ < 1 under the assumptions in force on V{9). 

Now in addition to j^j < 1 we choose arbitrarily a constant ct > and consider those z for which 
logd^l) < —{k — (T)log{n)/n (this is the interesting case, since according to Corollary 13.21 we are excluding 
a zero-free annulus near the unit circle whenever a > S where S is the arbitrary positive parameter in the 
statement of Corollary 13. 21 note that both S and a may be taken to be arbitrarily small). Since the function 
r"| log(r)|'"' achieves its maximum value when | log(r)| is proportional to l/n, we have 



(162) max r"|log(r) 

log(r)<-(fe-cr) log(Ti)/n 



(A:-a)*'n'^-^'=(log(n))* 



r"|log(r)|'= 

log{r} = -{k-a) log(n)/ri 

Since we are assuming here that A: > 2, ct > may be chosen small enough that both error terms in H161|) 



may be replaced by o(n' 



-(fc+i) 



as n ^ oo uniformly for log(|z|) < — (/c — a) \og{n)/n. 



It remains to evaluate the explicit integrals in (|161|l by Laplace's method. Letting a = fcora = fc + lwe 
consider 



(163) 



YaA^) ■■-- 



M\s\,0,){-log{\s\)r 



ds 



2-./2g.ej 



a;"/i(a;,6lj)(-log(a;))'^ 



dx . 



t/2 
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Expecting the dominant contribution to come from the neighborhood of a; = 1, we write 
(164) 



YaA^) 



'7i(l, 



1 — ze 



2-«/2 



a;"(-log(a;))"rfa; + e™''^- 



2-'/^ 



h{x, 9j) 



HI, 



1 — ze 



Finding a common denominator and extracting from Yaj{z) a factor of e*"^^(l 



x"(-log(a;))"da;. 
)^^, we see that 



^in9j 



YaA^) 



(165) 



hii,e,) 



t/2 



x'\-\og{x))''dx 



ze-''^^ {h{l, Oj) ~ h{x, Oj)) + h{x, 9j) ~ xh{l,9j) 



a;"(-log(a;))"da; 



If z does not approach the point e'^^ , then the fraction in the integrand of the second integral is easily seen 
to be bounded by a multiple of — log(a;) that is independent of n. More generally, if z is allowed to approach 
the point e'^^ as n — > cw, then we note that the integrand is analytic in x, and the path of integration may be 
deformed in such a way that throughout the path of integration \x — ze^*^^ | is bounded away from zero by a 
quantity that is proportional to log(n)/n because log(|z|) < —{k — a) log(n)/n, and therefore the fraction in 
the integrand of the second integral is bounded by a multiple of — log(a::) that is proportional to n/ log(n). 
Therefore, 



(166) 



YaA^) 



'^^h{l,0j) a\ 



l + O 



1 



1 — ze ^"^ 77,°+! \^ \^log(n)^ 

where the error is uniformly small for log(|z|) < —{k — cr) log(7T,)/n. Since EkV{l,9) = V{9), we have 
h{l,9) = e-^o+*f^(«). It follows that 

e 



(167) 



1 

n 



W^L.12('^',^') 



dA' 



jfc-ig-yo 



2j^j^fe+ 



■ „ . . ,, , „i(n+l) 

1 Z^ J pie J _ 



-(fc+i) 



2-^/2<r'<l 

holds uniformly in the region log(|z| 



< — (fc — (t) log(n)/n, and therefore by p38|) so does 



^(fc+i) 



(168) -4^) = ^"e----(V-)e-^-(^^«)B(log(.)A) + ^^^^^H^ , 

This concludes the proof of Theorem 13.51 

4. Exponentially Varying Weights 

4.1. Asymptotic formulae for 7r„(z) and 7„ in the varying weights case. In this section, we consider 
weights of the form 



-nV{e) 



for all 9eS^ 



(169) 

where y : S*^ ^ M is a given real-valued function of period 27r. The weight (|169|) varies exponentially 
according to a parameter tt., and for each n we may associate with the corresponding sequence of monic 
orthogonal polynomials 7ro(z), 7ri(z), 712(2:), . . . and normalization constants 70, 71,72, ... . Properly speaking, 
these quantities depend on the parameter n, and we should invent notation to express this dependence, such 
as TTj {z). We will not introduce this cumbersome notation. However, the reader should take note of this 
dependence. The limit of interest here is to study the behavior of the particular monic polynomial iTniz) 
of degree n in this system along with its normalization constant 7„, in the limit n ^ 00. Thus the large 
parameter n enters simultaneously into the degree of the polynomial and also into the weight, and we are 
studying the asymptotic behavior along the diagonal of a doubly indexed sequence. 

The asymptotic behavior in this limit is governed by the function V, along with some associated functions. 
First, recall the analytic function N{z) defined by H2U|I for \z\ > 1, which is associated with the negative 
frequency component of the Fourier series of V(9). Now define 



(170) 



k{9) 



N{e"')-N{e'0) 



ni9), 
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where n{d) is the periodic function defined in H21|) . Both functions k and fl are as smooth as V is. 
Our asymptotic results in this case are the following. 

Theorem 4.1. Let p > be a fixed integer. Suppose that (f>{6) = e^"^^^-* where V : S^ ^^ R is of class 
(^fc-ia^^i-j jg^ some k > 2p + 2. If k'{9) is strictly positive, then for each p > 1 there is a constant Kp^p > 
such that 



(171) 



sup 

\z\>p 



dP 



dzP 



7r„(z)z-"e-"^(^) 



1 



< K, 



log("-) 

P.P„2(fc-l) 



holds for all sufficiently large n. 

Theorem 4.2. Let p > be a fixed integer. Suppose that (j){9) = e^"^'^-* where V : S^ ^ S. is of class 
(^fc-ia^^i-j y^j, some k >2p + 2. If k'{9) is strictly positive, then there is a constant Kp > such that 



(172) 



sup 

\z\>l 



dP 



dzP 



r„(z)z-"e^"^(") - 1 



<K, 



P„fc-2p-l 



holds for all sufficiently large n. 

< Remark: As in the fixed weights case (see (|42|l '). a special case of the estimate (|172() is the following 
estimate (holding under the same conditions) 



(173) 



sup 

-7r<9<7r 



-le 



dO 



I i0\ -inS -nNie'") 

7r„(e )e e ^ ' 



<K 



. log(n) 



-.k-2p-l ' 



asymptotically characterizing 7r„(z) on the unit circle. > 

Once again, the uniform nature of the convergence on the circle allows us to prove the following mean 
result (the proof is the same as for the analogous result in the fixed weight case) . 

Corollary 4.1. Let p > be a fixed integer. Suppose that 4>{0) = e^"^*-^' where V : S^ —> M. is of class 
(^fe-i,i(-^i^ wit/i fc > 2p + 2. If k'{0) is strictly positive, then 



(174) 



hml.t^ 



n — ^oo jiP 



Tn||0 



The next result concerns the asymptotic behavior of the polynomial 7r„(z) for z inside a closed annular 
region whose outer boundary is the unit circle. 

Theorem 4.3. Suppose that 4>{e) = e-"^^''' where V : S^ -^ M. is of class C'^-^'^iS^) for some k > 2. 
If K,'{9) is strictly positive, then for each p satisfying < p < 1 there are constants K^ > such that the 
estimates 

(175) 

and 

(176) 



sup 

P<|2|<1 



7r„(z)e 



-nN{l/z) 



£n^inEk^{rM) 



< K 



- log(n) 



sup 

\z\<p 



7r„(z)e 



-n.Af(l/z) 



< 



^ 

'n"' ^ 



hold for all n sufficiently large. 



An immediate corollary is that there exists an annulus inside the unit circle that asymptotically contains 
no zeros. 

Corollary 4.2 (Zero-free regions). Suppose that <p{9) — e~"^^^^ where V : S^ ^ R is of class C''~^'^{S^) 
for some k > 2, and suppose that k'{9) is strictly positive. Let 6 > be an arbitrarily small number. Then 
there are no zeros o/7r„(z) in the region 



(177) U 

as long as n is sufficiently large. 



log(r) > - 



k — 1 — S\ log(?i) 



l + n'{9)J n 
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<i Remark: In the case of a fixed weight, Corollarv l3 . 2l established the existence of a zero-free annulus with 
outer radius 1, and whose inner radius depended explicitly on k, the degree of smoothness. We subsequently 
considered a family of weights for which wc could compute explicitly the behavior of the zeros, and showed 
that Corollary 13. 21 is sharp. Indeed, for explicit families of weights with k degrees of smoothness, some zeros 
achieve a distance of fclog(n)/n + log(log(n))/n from the unit circle. In addition, for these examples, a 
majority of the zeros approach a circle whose radius is 1 — (fc + 1) log(n)/n. For the case of varying weights, 
it is our belief that Corollary 14. 21 is similarly sharp. Moreover, it is to be expected that for some canonical 
family of examples constructed such that i^^^' has jump discontinuities, it should be similarly possible to 
obtain a very detailed asymptotic description of the zeros. i> 

<i Remark: The zero-free region is determined by observing that for there to be a zero of 7r„(z), 
^n^inEkCt(rfi) ^cm^f^ ]^q roughly the same size as log(n)/n'^~^. This is clearly not true for z on the unit 
circle, and investigating the size of 2:"e™^''"^^'"'^^ relative to log(n)/n''~^ yields the result. For \z\ < 1, 
^n^inEkn(rfi) jg gj-^aji not only because z" is small, but also because (as will be clear in S 14.21 below) we 
constructed g™Ekiiir,e) g^ |-]-^g^|- j-^j, 7. < i but 1 — r sufficiently small, 

(178) ^inE^nirfi) ^ Q /nn'(e)\ 

Thus our extension of the function fl plays a role in determining the zero-free region near the unit circle, o 
<i Remark: The quantity l + il'{6) appearing in Corollarv l4.2l is strictly positive because l + i^'{9) — k'{9), 
and k'{0) is strictly positive. It is interesting to note that this condition also guarantees that there are no 
gaps in the support of the equilibrium measure (see Appendix IX|) . The occurrence of a gap in the support of 
the equilibrium measure is heralded by the development of a zero of the function k' [6). Although Corollary 
14.21 would not apply if k!{9) vanished at some Oq, an intuitive consideration of the set defined in H177|) 
indicates that near ^Oi the zeros are pushed away, further from the unit circle. As a gap develops then, one 
might expect the zeros to accumulate on a contour approaching the unit circle, but with a gap aligned with 
the gap in the support of the equilibrium measure. We will not carry out such an analysis here, but clearly 
the methods outlined here can be adapted in this direction. \> 

Just as Theorem 13.31 leads to Corollary 13.31 (|176|l from Theorem 14.31 yields the following result. (Recall 
from the definition l|2()|l that N{z) ^ as z — > 00, and so iV(l/z) ^ as z -^ 0.) 

Corollary 4.3 (Varying recurrence coefficients). Suppose that (f){9) ~ e~"^'^) where V : S^ ^ M. is of class 
C'^~^'^{S^) with k > 2. If k'{9) is strictly positive, then there is a constant K > such that the estimate 

(179) |a„| < -^ 

holds for sufficiently large n. 

Finally, we have the following result concerning the asymptotic behavior of the normalization constant 
7„ > 0, defined such that ||7n7r„(z)||0 = 1, with 0(6*) given in the varying weights case by p69|l . 

Theorem 4.4. Suppose that V is a real function of class C^^^'^{S^) for some k > 2. If k'{9) is strictly 
positive, then there is a constant K > such that 

holds for all n sufficiently large. 

< Remark: At first glance, the asymptotic formulae H18U|I and H179|) may appear to be inconsistent with 
the identity l|33() . However, in H180|l . 7„_i is the norming coefficient of the {n — l)-st degree polynomial 
orthogonal with respect to the n-dcpendcnt weight e~"^. To verify 1)33(1 in the varying weights case, one 
would have to compute the asymptotics for 7„ as well, rather than merely replacing ri by n + 1 in p80() . > 

4.2. The d steepest descent method for exponentially varying weights. One of the main points 
of this paper is that while the limit appropriate for exponentially varying weights of the form 1(1 69() lies 
for the most part beyond the reach of classical techniques applicable for fixed weights (meaning primarily 
the approximation of (f){9)~^ by positive trigonometric polynomials), analysis of this limit by means of the 
d steepest descent method presents almost no further difficulty beyond the analysis carried out for fixed 
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weights in § 121 To illustrate the ease with which many of the techniques carry over from the fixed weights 
case, we now outline the analogous calculations for the varying weight (|169|l . 

As in the case of fixed weights, wc begin with the matrix M"(z) solving Ricmann-Hilbcrt Problem 12.11 
and introduce a sequence of explicit transformations. In order that Riemann-Hilbert Problem 12.11 indeed 
characterizes the polynomial of degree n orthogonal to all lower degree polynomials with respect to the 
weight H169|) . we make the following assumption. 

Assumption 4.1. V is a real continuous function on the circle that, for some exponent v G (0,1] and for 
some constant K > 0, satisfies a uniform Holder continuity condition \V{92) — V{9i)\ < A'|02 ^ ^iT- 

4.2.1. Conversion to an equivalent d problem. Solution of the d problem in terms of integral equations. Now, 
define S"(2:) in terms of M"(z) as follows: 



pnVb<T3/2]yjTi/^\g-riVo<T3/2 



^ — n—nN{z) 











(181) S"(z) := < 



^nVoa-3 /2-^n / ^\^-nVoa-3 /2 



for \z\ > 1 

for \z\ < 1 . 

Since N{z) — > as z — > cxd, it is easy to see that as z ^ oo, S"(z) — » I. Moreover, S"(z) is analytic 
for \z\ ^ 1, and its boundary values S"(z) (respectively S"(z)) taken on the circle S^ from the inside 
(respectively outside) satisfy: 

„mre(e) 



g-«W(l/z) 





p"JV(l/^) 




(182) S!^(e*'') = S'!(e*'') 

e 

Note that according to Assumption 14. l| k(6) also satisfies a uniform Holder continuity condition with expo- 
nent v. 

As in § 13.2.11 we have available an algebraic factorization of this jump condition: 

/ 1 o\/oiW 1 

(183) S!^(e*^) = S:!(e*^) 

\ £""""(*') 1 / \ -1 / \ e™"'*) 1 

To take advantage of this factorization we need to extend the functions e^™''^^) from S^ to an annulus 
containing S^. To have continuity of these extensions, wc make the following assumption. 

Assumption 4.2. The function V is of class C'^^^{S^) for some fc = 1, 2, 3, . . . . 

With this assumption, which provides new information only if fc > 1, we can extend e^™'*^^' as follows. 



Writing e 



±inK( 



^±m0^±mn(e)^ we extend the factor 6=*=™^ analytically as 2=*=". On the other hand. 



r2 : S^ ^ M is a well-defined function on the circle, and therefore we may apply the extension operator Em 
defined in ij^l to fl, resulting in an extension of e^'"'^^W to the domain R^ \ {0} as a continuous function 
g±m_E„n(r,e) fpj. g^j^y j^ ^^ ^Yie range 1 < m < k. 

Let e > be the radius parameter of the annular domains A± (see Figure|Hll. Recall the "bump" function 
B defined in S 11.51 We define a new matrix unknown TJJ^ ^(r, 9) by setting 

1 



(184) T;,Jr,e) 



S"(z) 

S"(z) 
S"(z), 



B(log(r)/e)r^"e"™^e-™-^'"'^^('''^) 1 

1 

-B(log(r)/e)r"e"''e™^'""('''^' 1 
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for z = re'^ G A. 



for z = re'^ £ A 



for z = re'^ ^ A+ U A. 



Note that the presence of the factor B(log(r)/e) ensures that T'!^^{r,9) may be continuously extended to 
the outer boundary of A_ and the inner boundary of ^4+; that is, T^ ei'f'yf^) is a continuous function for 
r ^1. 

< Remark: Here our analysis rests upon extending S1(0) from the circle into domains A±. By contrast in 
our treatment of the fixed weights case in S 13.21 it was more convenient to extend the function V{9), which 
is related to il{9) by a Cauchy transform. t> 

li m ~ 1, then it is easy to see that EmQ,{r,9) ~ Eiil{r,9) = n{9) for all r > 0, so the off-diagonal 
matrix elements in S"(re*^)~^T"^(r, 0) are bounded in magnitude by e~"l'°s('")l. Therefore, in this case 
()184|l represents an exponentially near-identity transformation when n is large and | log(r)| is not too small. 
We want to ensure that a similar situation prevails when m > 2 as well. Since 1 < ttj < fc, the consideration 
of m > 2 requires that fc > 2. In this case Assumption 14.21 implies that k{9) is continuously differentiable, 
and hence so is il{9). The crucial conditions for H184|l to be a near-identity transformation when n is large 
are that Q'{9) + 1 is strictly positive and that the parameter e is chosen small enough, as the following lemma 
shows. 

Lemma 4.1. Suppose that n(9) is a real function of class C™^^(S'^) for some m>2, such that Q,'{9) + 1 
is strictly positive. Then there exist constants eg > and fi > such that whenever < e < eoi 

(185) g-»B™n(r,e) < ^i-p ^ j^^ ^11 9 and forl<r<2', 

and 

(186) g«B™n(r,e) ^ ^^-i ^ j^^ ^^ q ^^^ j^^ 2"^ < r < 1 . 



Proof. By hypothesis, we have 
(187) 



inf Vt'(9) = r - 1 , where r > . 



Note that log |eT'-E™f^('-.e)| ^ ±3(£;,„rj(r, 0)). From ©, we have 
(188) 

^E,M{r,9)) = J2 (-l)"-^' (2 +1), l"g(0'"+^ = -log(r) 



p=0 



p=l 



{2p+\)\ 



where P{m) = {m — 2)/2 if m is even and P[m) = {m — 3)/2 if m is odd. Since fl has m — 1 continuous 
derivatives we have 

Pirn) sup |17(2P+1)((?)| 



(189) 






p=i 



for all r satisfying | log(r)| < elog(2) where e > 0. Therefore, eo > may be chosen sufficiently small that 
the inequality < e < eo implies that 

P(r 

(190) n'{9)+J'(-^^p- 

p=l 

holds whenever | log(r)| < elog(2). Then, we have 



Q{2p+l)('fl^ J- 

n\9) + y: i-ir L , A? logM^^ > J - 1 



(2p+l)! 



(191) 

and 
(192) 



log 



log 



-i_E,„0(r,e) 



JE,„n{r,0) 



< - log(r) 



< log(r) 



I-l 

2 



holding for 1 < r < 2"^ 



holding for 2 ^ < r < 1 



The estimates (|185ll and (I186II thus both hold with the choice fi :~ r/2 > 0. 

Thus, we are led to propose the following assumption. 
Assumption 4.3. If V is of class C^{S^), then k'{9) defined by p7U|) is strictly positive. 
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D 



This assumption is not an explicit assumption on V{0), however it is possible to give conditions on V{9) 
that are sufficient to make Assumption 14 . 31 hold . For example, since for V G C^(5^), we have 

oo 

(193) K'(e) = 1 - 2;^ j|^,| cosO'^ + arg(VS)) , 

where we recall the Fourier coefficients of V defined by (|19() , the condition 



(194) 



E I 



\V,\<1 



is sufficient to guarantee that Assumption 14.31 holds. If in fact V G (7^(5^), then the convexity condition 
(195) 






also guarantees that k'{6) is strictly positive. Our proof that the condition H195|) implies that k!{6) is 
strictly positive makes use of certain aspects of logarithmic potential theory and is given in Appendix IXI In 
Appendix^ we show that neither are the sufficient conditions (|194|l and H195|) equivalent, nor does either 
condition imply the other. Rather, the two conditions are independent and thus complement each other. 

<\ Remark: In the approach taken in |2| , where a particular case of a varying weight in which the function 
V{9), and hence k,{9), is an analytic function, the analytic extension Eoo^{r,0) is used, and the "bump" 
function factor S(log(r)/e) is omitted. The latter omission has the effect of introducing exponentially small 
jump discontinuities in T", ^(r, 9) across the circles E±. o 

Since for m = 1, . . . , fc we have Em^{^, 9) = 0(0), the matrix TJ^(r, 9) satisfies the jump condition 



(196) 



limT:^,,,(r,0)=limT:;^,(r,( 



1 



-1 



To remove this jump discontinuity, introduce one further change of variables: 

r T" (r,0), forr>l 



(197) 



37aAr,e) 



T" (r,0) 







1 



for r < 1 , 



At this point, we can relate J^_e(r, 9) directly to M"(z). Combining ifTM]) . ifTMj) . and (fW7| . we have by 
definition 



(198) j;;jr,( 



^nVoa^ /2y^n /^\^~n.Voa3 /2 







,-nN{l/^) 



\ 



^nVaa3/2y^n/^\-nVQa3/2 



for < 7' < 1 

griJV(z)^^Q„/j,Wg\g-inB„0(r,e) ^nN (z) ^n ^inO 

for r > 1 . 

Given the assumptions in force, this matrix function is clearly continuous throughout the plane. To determine 
its deviation from being an analytic function in the regions r < 1 and ?' > 1, we need to control a derivative, 
and consequently we make the following assumption. 

Assumption 4.4. The Junction V is of class C''~^'^{S^) for some k > 1. That is, V'-''~''^\9) is Lipschitz 
continuous. 
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Then, for 1 < m < k, we see by direct calculation that 



gnVo(T3/2-|y|-n/^-ig-nV'o(T3/2 



(199) dJ^^Jr, 







giiVo cr3/2-|y|-n /^-.g-nVbo-a 



/2 



e"^(i/^V"e™*5 [B(log(r)/e)e™-^""(''''')] 

for < 7' < 1 


for r > 1 



Eliminating M"(2) in favor of JJ^ ^(r, 6') using (|198|l again yields 

(200) 9j;^_,(r, 61) = J;;,,,(r, 6i)X;;,_,(r, 0) , for r 7^ 1 and almost all 6 e S^ 

where 

r"e™^9 [B(log(r)/e)e"-^-"(''''')] 



(201) 



^mA^ 











m,e\ ' 











for < r < 1 



for r > 1 



^-n^-ineg [5(iog(r)/e)e-™-^"^^(''^^)] 

Then it is easy to see that JJJj ^{r, 9) satisfies the following d problem. 

d Problem 4.1. Find a 2 x 2 matrix \J{r,6) with the properties: 

Smoothness. XJ{r,9) is a Lipschitz continuous function throughout M.'^ . 
Deviation From Analyticity. The relation 

(202) aU(r,0)=U(r,0)X:;,,(r,0) 

holds for all points in R^ with the exception of a set of measure zero. The matrix XJ^ ^.(r, 9) is defined 
almost everywhere by (|201|l and is essentially compactly supported. 
Normalization. The matrix \J{r, 9) is normalized at r = 00 as follows: 

(203) lim U(r, 61) = I . 

r — ^oo 

In exactly the same way that Proposition 13 . II was proved, we have the following. 



Proposition 4.1. Suppose that 1 



- „-«v 



where V : S^ 



is of class C-^-^iS^) for some k = 1,2,3,. 



Then for all n = 0, 1, 2, 3, ... , for m = 1,2,...,/;:, and for all e > 0, the matrix XJ^ ^.(r, 0) is well-defined 
almost everywhere by H2U1|) and d Problem \4-. l\ has a unique solution, namely U(r, 0) = JJ^ ^(r,9). 

Also, the proof of Proposition 13 . 21 carries over to the context of the varying weight p69|) in the following 
form. 



Proposition 4.2. Suppose that 1 



- „-«v 



where V : S^ 



is of class C-^-^iS^) for some k^ 1,2,3,. 



Then for all n = 0, 1, 2, 3, ... , for m = 1,2,...,/;:, and for all e > 0, the matrix X"^ ^(r, 9) is well-defined 
almost everywhere by (|2U1|) and the corresponding solution U(r, 0) = 3'^^{r,9) of d Problem \4-l\ satisfies 
the integral equation 

U(r',0')X" (r',0') 



(204) 



U(r,0)=I-iy|. 



■dA' 



z — z 

where z — re , z' = r'e , and dA' is a positive area element dA' = r' dr' d9' . The integral is taken over 
the entire plane. 
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4.2.2. Asymptotic solution of the integral equation. Estimates of 3"^ ^{r,9) and its derivatives for large n. 
As in § 13.2.21 it is possible to characterize JJ^ ^ (r, 6) by analyzing the integral equation (|204|l as long as n is 
large enough. However, in the context of the varying weight H169(l we will require the monotonieity condition 
expressed in Assumption ^31 ^^^d that the radius parameter e be taken sufficiently small for each admissable 
given V that the conclusion of Lemma [4. II holds. Moreover, the exponential character of the varying weight 
()169|l suggests that, by comparison with the analysis presented in § 13.2.21 this approach is only fruitful in 
giving the same degree of control on JJJ^ ^(r, 9) as was achieved for H^ ^{r, 9) if V{9) has more smoothness 
than was required in Proposition 13 . 41 In particular, we will require that the following assumption holds. 

Assumption 4.5. The function V is of class C''~^'-^{S^) for some k > 2. 

In other words, to achieve the same convergence rates as in the fixed weight case, one more derivative of 
V will be required in the exponentially varying weight case. 

As in i) 13.2.21 we need some bounds for XJJj g{r,9) and its derivatives. 

Proposition 4.3. Suppose that V : S^ ^ M. is a real function of class C''^^'^{S^) for some k > 1 for 
which fl'{9) + I is strictly positive, that m is an integer satisfying 1 < m < k, and that e > is sufficiently 
small. Let the integer D he defined as D := min(fc — to, m, — 1). Then, the matrix function XJ^ ^ is of class 
C^"^'^(R2 \ {0}) if D > 0, and of class L^{R'^ \ {0}) if D = 0. Moreover, if a and (i are nonnegative 
integers such that a + f3 < D, then there are constants /i > and Cm,e > such that for all n the estimate 

Qot+p 



(205) 



-X" ir 9) 



< C^/^ni+'^e""^! '°eMI I log(r)|™-i-" ^ nP\ log(r)|? 

p=0 



dr°'d9l^' 
holds throughout the region |log(r)| < elog(2) containing the essential support of "K.^ ^[r,9). 



Proof. This proposition is proved in almost the same way as Proposition l3.3l In this case, it is essential to 
recall Lemma [4.11 which provides the constant /i > and thus the exponential decay of the term e~"^l iog(i-)l 
for e > sufficiently small. Also, the initial application of the d operator in the definition H201(l yields a 
factor of n that does not appear in the proof of Proposition l3.3l D 

<i Remark: This result should be compared with ProDOsition l3.3l The only important difference between 
the estimate (|2Q5|I and (|96|l is the presence of an additional factor of n. > 

The proof of Proposition l3.4l with references to Proposition 13.31 replaced by references to Proposition l4.3l 
applies to the asymptotic estimation of 3^ ^{r,9), with the only important difference being an additional 
factor of n. This results in the following. 

Proposition 4.4. Suppose that (f) = e~"^ where V : S^ ^tM. is of class C^~^'^{S^) for some fc = 2, 3, 4, . . . . 
Let the integer m lie in the range 2 < m < k and fix e > sufficiently small. Define the integer D :~ 
min(fc ~m,m — 2) > 0. Then, for all n > the matrix XJJ^ ^(r, 9) is well-defined almost everywhere by (|2U1|I . 
and for all n sufficiently large, 3"^ ^(r, 9) is given by a Neumann series 

(206) 3Z,Ar,0)^l+ iX:iJ) {r,9) + ( A-^^,, o x:,\J) (r, 0) + • • ■ 

which converges in the norm ||| • |||^, where the double integral operator X"^ ^ is defined by 

1 f rF(r',9')X.''Jr',9') , 

(207) {XZ,,F){r,9):=--J J ^ ^, ^-^^ ' dA' . 

In particular, if D ~ then JJ^ ^ lies in the space L°°(M^), and if D > then J^ ^ lies in the space 

C^~^'^(]R^) and IHJmelllf) is finite. For all integer p in the range < p < D, the following estimates hold 
for sufficiently large n: 

(208) \\\r^,e-mp<C^ll^^, 
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where Cm!e is a positive constant. Furthermore, for each p > 1'^ and j or all integer "p in the range < p < D, 
the following estimates hold for sufficiently large n: 

Qa+I3 



(210) 



E 

a+l3<p 



sup 

-ir<6l<7r 
|log(r)|>log(p) 



dx^'dyP 



[Jr„,,(r,0)-I] 



< &!). 



1 



™,P^m-p-l ' 



(211) 



E 

a+/3<p 



sup 

-7r<e<7r 
log(r)|>log(p) 



Qa+13 



dx^'dyP 



[r^^,{r,e)-i-{x::,Mr.e)\ 






log(?t) 



Mp) 



where Cm'p is a positive constant. 



< Remark: Note that in H2Q8() . due to the extra factor of n introduced into the estimates by Proposi- 
tion ^31 convergence of the Neumann series in the ||| • |||p norm foUows provided that m — p — 1 > 0, and 
so wc inay only consider m > 2. In order to uniformly control p derivatives of JJ^ ^{r, 9), Proposition 14.41 
requires that m should lie in the range 2+p<m<k~p, and therefore in order for there to exist suitable 
values oi m, V : S^ ^ R should be of class C'^~^'^(S'^) for some fc > 2p + 2. > 

< Remark: Assumption 14.31 that V{9) is such that k{9) is a strictly increasing function of the angle 
9, ties together two crucial aspects of our analysis. First of all, as can be seen from the matrix factors 
involved in the change of variables (|184|) between S'"{z) and TJ^ ^(r, 9), the inequality k'{9) > is precisely 
what makes this a near-identity change of variables in the regions f2±. On the other hand, it would not 
suffice to replace k'{9) by another unrelated positive quantity, because we require of whatever extension we 
choose some degree of vanishing of the d derivative at \z\ = 1 in order to control the d problem (that is, to 
sufficiently bound the operator X^ ^). This vanishing is built into the extension operators we have defined in 
(|5|) by the key property Q. Moreover, any smooth extension of k{9) that has a vanishing d derivative on the 
unit circle will behave similarly near the unit circle. So wc conclude that while the analyticity of k(9) is not 
important, the monotonicity of this function is crucial. To make an analogy with the asymptotic analysis 
of oscillatory exponential integrals, the d steepest descent method resembles Kelvin's method of stationary 
phase more than it does the saddle point method. On the other hand, it is important to note that the d 
steepest descent method remains fundamentally a method of deformation into the complex plane, and is not 
based on integration by parts^, the Riemann-Lebesgue Lemina, or related arguments of harmonic analysis. 
For an approach based on the latter, see Varzugin [5^1. o 

4.3. Proofs of theorems stated in § 14.11 The proofs are generally based on expressing 7r„(z) = Af[\(z) 

in terms of explicit functions and JJJj ^(r, 9) by means of 1198|l . Then, one applies Proposition 14 . 41 to control 
J^ g(r, 9) —I and its derivatives. 

4.3.1. Asymptotic behavior of iTn^z) for \z\ > 1. Proof of Theorems \4.1\ and \4.'^ From (|198|) wc have the 
following exact representation for 7r„(z) valid for \z\ > 1: 



(212) 



T^niz) 



^nN{z) 



\z\>l. 



^"■/™,e.ii(^,^) - B{\og{r)/e)e-"^^-^'^^^'^Jl^^^,^{r^ 

Here e should be taken to be sufficiently small. It follows that 

(213) ^„(z)z-"e-"^(^) - 1 = [Jl^.siir. 9) - l] - Bilog{r)/e)z-e-^"^-''^^''^ J:^^,^,,{r, 9), \z\ > 1 . 

If jO > 1 is fixed, then perhaps by making e > smaller yet, it can be arranged that _B(log(r)/e) = whenever 
\z\ > p. In this case, we have 



(214) 



dP 
dzP 



7r„(z)z-"e-"^(^) 



d''[j:^.,,ni^,0)-l] 



\z\>p>l. 



Properly speaking, we do not rely on integration by parts (or more generally, Stokes' Theorem for the d operator in the 
plane — in as much as this can be considered a generalization of the standard d inversion formula) to establish the existence 
of an asymptotic expansion. However, such methods are useful in the detailed analysis of individual terms in the expansion. 
This technique was used, for example, in the proof of Theorem 13. 51 
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Using the estimate (|211f) from Proposition 14.41 in the case m = k — p, and the fact that X", ^ (r, 0) is 
off-diagonal, we sec that for some constant C > 0, 



(215) 



sup \dP [4%.Mi('^. ^) - 1] I < G^^^ 

r>p>l 



Now since JJ!_„ ^ ui^) — 1 is analytic for \z\ > p, and tends to zero like z~^ a.s z —^ oo, Cauchy's Theorem 
for an exterior domain, together with H215|l for p = 0, proves (|171|l . 

To prove H 173(1 . we fix e sufficiently small and consider 1 < r < 2^/^ in which case i3(log(r)/e) = 1 so that 



(216) 



dzP 



Kr < 2'/2 . 



7r„(z)z-"e-"^W -l\^dP [J;^,,.n(^e) - l] - d^ [z-"e-'^"''-''^^^''>j:^,,^,,{r,9) 
Again, since XJJj ^(r, 9) is an off-diagonal matrix, we see from (|209|l that there is a constant C > such that 

log(n)2 



(217) 



sup|9P[j,%,,^n(^,^)-l]|<C'- 



,2(fc-2p-l) 



As in the fixed weights case, the dominant contribution comes from the remaining terms on the right-hand 
side of H216|l . Using H2()8(l from Proposition 14.41 to see that v?d^~^ J^^ i2(''!^) is of order log(n)/r 



m—p—l 



and taking the best case of tti = fc — p, we then find that 

dP 



(218) 



sup 

-7r<e<ir 
l<r<2=/2 



dzP 



■Kn{z)z 



-n —nN{z) 



<K 



, log(n) 



-,fe-2p-l 



This proves (|173|l . Now the maximum modulus principle applied to 7r„(z)z "e "^(^) implies (I172|l . 

4.3.2. Asymptotic behavior o/7r„(z) for \z\ < 1 and o/7„_i. Proof of Theorem \4.S\ and Theorem, \4-.4\ The 
proof of Theorem 14. 41 is based on the identity 7^„i = — M^;^(0). Using (|198|l . we see that 



(219) 



7Lie""''° = ^™,.,22(0,e): 



and thus Theorem 14.41 is proved by applying the estimate (|211|) from Proposition 14.41 in the case of p = 
and m = k. 

Proving Theorem 14. 31 begins with the exact formula 



(220) 



7r„(z) = -e"^(i/^) j;\^,.i2(r, 9), < r < 2" 



which follows from ((TOS|l using 7r„(z) = Af{\(z) and the fact that S(log(r)/e) = for r < 2-^ Using l(?Tn|) 
from Proposition 14 . 41 with p = and m = k then completes the proof of (|176() . 
For the proof of (|175|) , one begins with the following formula for 7r„ , 



(221) 



^„(z) = -e"^(i/^)j;;.,,.i2(r,0) +e"^(i/^)e-™^^"(^'«)j;^,,, ,i(r,0). 



which again follows from (|198|l . and the fact that i?(log(r)/e) = 1 for r > 2 ^/^. Using (|208|l from Propo- 
sition 14.41 with p = and m = k then completes the proof of 1)175(1 , and this completes the proof of 
Theorem O 
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Appendix A. Logarithmic Potential Theory of Orthogonal Polynomials on the Unit 

Circle 

A more general and systematic strategy for extracting asymptotics of M"(z) in the exponentially varying 
weight case when 0(6*) is of the form (|169|) is the following. First introduce a function g{z) (to be determined) 
satisfying g{z) ~ log(2;) as z — > oo and that e^^^' is analytic for \z\ ^ 1, taking continuous boundary values 
on the unit circle E. Then one converts Ricmann-Hilbcrt Problem 12.11 for ]V["(z) into one with identity 
asymptotics as — > oo by the change of variables 

(222) M"(z)=Y"(z)e"s(^)'^3^ 

resulting in a new unknown matrix Y"(2). In addition to the normalization condition lim^^oo Y"(z) = I, 
Y"(z) satisfies the jump condition 

(223) Y!^(e''')=Y!!(e^'') , ior 6 € S\ 

where the subscript "+" indicates the boundary value taken from within the circle, and "— " indicates the 
boundary value taken from outside. Without any loss of generality, we may represent g{z) in the form of a 
complex logarithmic potential 

(224) g{z)= f Lg,{z)4,{e')de' 

J —IT 

where for each 6 £ S^ wc consider the function Lg{z) := log(z — e'^) to be real for z — e*^ sufficiently large 
and positive real, with the branch cut from the point z = e** in the clockwise direction along the unit circle 
to the negative real axis, and then along the negative real axis to z = —00. We also require that 

(225) / ilj{e')de' = l 



in order to satisfy the required normalization condition g{z) ~ log(z) as z — > 00. 

Additional conditions may now be placed on g, or equivalently on ip, in order to make the jump condition 
(|223|l for Y"(z) asymptotically tractable. One key condition is that the function ip should be real- valued. 
From this condition it follows that g+^z) + g-{z) — log(z) has a constant imaginary part for \z\ = 1, as the 
following argument shows. When \z\ = e'^, the identity d/d6 = izd/dz yields 

(226) -^^{g+iz)+g_iz)-\ogiz))^^{zg'4z) + zg'_iz)-l) . 

Using H225|l and differentiating 12241) under the integral with respect to z for \z\ ^ 1, we obtain 



(227) ^cj(g^(^) + g_(^)_log(z)) = -limK 

dd ej.0 






^ (e«e' _ eze)2 _ ^2^ 
Assuming reality of "0 we can pass the real part under the integral and thus arrive at 

/ooo^ ^c-/ /N / ^ , / ^^ ,• T e^ + 26^ + 2e2 cos(C) - 46^ cos2(C) ,,n:^,a: 

(228) -5 (.9,(^) + 5-(^) - log(.)) = - hm|^ 4 + .^ + (4.^-8)cos(C) + (4-4.^)cos^(C) ^'^'^ '' 

where C = 6' — 6. The numerator of the fraction in the integrand is clearly uniformly bounded by a quantity of 
order e^, and the minimum value of the denominator is achieved when cos(C) = (2 — e^)/(2 — 2e^) G (—1, 1) 
yielding a minimum value that has the asymptotic expansion e^ — e* + O(e^) as e [ 0. Consequently 
the fraction is uniformly bounded independent of e. Moreover, the fraction is easily seen to tend to zero 
pointwise in the limit e j for ^ 7^ 0. A dominated convergence argument therefore shows that the limit on 
the right-hand side of (|228(l is zero. Sec also Pj, where it is also shown that the constant value is exactly 
9(5+ (z) -|- 5-(z) — log(z)) = TT. Another consequence of assuming that 1/) is real comes from noting that 

(229) 9+{e'^) - 9-{e'% = 27ri / ^p{0')de' . 

Je 

It therefore follows that g+{z) — g-{z) is purely imaginary for \z\ — 1. 
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Since 
(230) 5R(<?+(e'^) + g-ie''^)) = 2 / log e'^ - e'^' tPid') d9' , 



one is led to seek ip so that the circle is split into intervals of two different types: 
Bands: For 6 in a band /, il]{0) is real and positive, and 

(231) 2 I log e'^ - e^' ^{B') dO' ~ V{e) =£, for 6* in a band , 

where £ is a real constant (the same constant for all bands — in fact if there is only one band it 
turns out that i = — Vo). Thus, in a band the jump condition H223|l for Y"(z) takes the form 

(gin(K/(0)-7r) ^n{i+iTT) \ 

j , for e in a band / , 

where ki{9) is a strictly increasing real function. 

Gaps: For in a gap F, ip{0) = and we have the strict inequality 

(233) 2 / log e''^ - e'"' ipiO') d9' - V{e) < £ , for 61 in a gap . 

J —IT 

Thus, in a gap, the jump condition (|223|l for Y"'{z) takes the form 
/ g«i(Kr-7r) ^nie+tTT) . exponentially smaU \ 

(234) Y!f (z) := Y!!(z) , for 9 in a gap F , 

\ g-m(Kr-7r) / 

where Kr is a real constant. 

The alternative conditions (|231() and (|233|l are exactly the Euler-Lagrange conditions for the minimization 
of the weighted logarithmic energy 

(235) E[ij]:=l I iog(—l-^\^{e'Uie)de'de + l v{e')iie')de' 

over all probability measures il){9)d0 > supported on the unit circle. The constant i is the Lagrange 
nuiltiplier introduced to enforce the constraint (|225|l . The minimizing measure with density V'(^) is called 
the equilibrium measure. The connection of this extremal problem with the jump condition H223() through the 
Euler-Lagrange variational conditions suggests that logarithmic potential theory plays an important role in 
the asymptotic theory of orthogonal polynomials for general exponentially varying weights of the form (|169|l . 
The use of equilibrium measures for the asymptotic analysis of the Riemann-Hilbert problem associated to 
orthogonal polynomials on M with analytic exponentially varying weights was carried out in |1U[ [TT] . For 
orthogonal polynomials on the unit circle with analytic exponentially varying weights of a specific form, this 
was done in Pj. 

The transformation 

(236) Z"(z) = <^ 

[^ g-n£^3/2Y«(^)e"^'^3/2^ \z\>l 

leads to jump conditions of the form 

j , for z = e*^ in a band / , 

which is of exactly the same form as p82f) . and 

giriKr exponentially small 
(238) Zl{z) ^ Z'l(z) ( I , for z = e''' in a gap F . 
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Clearly, we also have the normalization condition Z"'(z) — > I as z — > oo. 

Now, if the function V{9) is such that the whole unit circle consists of a single band, or equivalently there 
are no gaps in the support of the equilibrium measure, then in fact Z"(z) — S"{z) as defined in 5 14.2. II and 
the analysis proceeds as in the main body of this paper. However, in the more general context — when it is 
only true that the support of the equilibrium measure consists of a finite number of disjoint intervals on the 
unit circle — a modification of the d steepest descent method described in § 01 is required. A more general 
method may still be based on the algebraic factorization H183|l of the jump matrix in each band; however 
the annuli A± must be replaced with a system of lens-shaped regions Aj_ adjacent to each band /. Since the 
variational inequalities become less effective near the band edges, a local analysis must be supplied to control 
the error. When such a method is developed, the conditions 1)194(1 or 1)195(1 on the function V{d) can be 
dropped as long as it is known that there are only a finite number of gaps in the support of the equilibrium 
measure. In |20| it is shown how to carry out such a program in the context of nonanalytic exponentially 
varying weights on the real line, where a convexity condition is known to guarantee the existence of a single 
isolated band [a,/3] C M. 

Let us show how the logarithmic potential theory described briefly above leads to the condition 1(1 9 5(1 
guaranteeing that the support of the equilibrium measure is the entire unit circle. If V is of class (7^(5*^), 
then for 9 in any gap in the support of the equilibrium measure the real function 

(239) ^{0) := 2 / log e''^ - e'"' ip{e') dO' - ¥{9) 



is also twice differentiable. In fact, a calculation shows that 

(240) ^"(9) - -2 /" W)d9' _ ^„ 

Note that this is not a singular integral since e*^ is assumed to lie outside the support of the equilibrium 
measure. Let us assume both the condition ((195(1 . and also the existence of a gap in the support of the 
equilibrium measure; wc will then derive a contradiction. Now, $(0) is continuous for — tt < 9 < tt. 
Therefore, since $(0) = ^ at both endpoints of the gap (according to 1(231(1 '). and $(0) < £ strictly in 
the interior of the gap (according to ((233(1 '). there must be a point 9 in the gap at which $"(0) > 0. However, 

since e** — e*^ < 2 for all angles 9' , we see that 

(241) $"(6*) < -J - V"{9) 

which is negative in view of the assumption l(195() , thus establishing the desired contradiction. 

Appendix B. Comparison of the Conditions (11941) and ((195() 

In this appendix, we illustrate by concrete examples how the two conditions ((194() and ((195(1 . while both 
sufficient for the prevention of gaps in the support of the equilibrium measure, are completely independent. 
Thus neither condition implies the other. 

First, consider the example V{9) — Acos{k9) where fc = 1, 2, 3, . . . and A G M arc parameters. A direct 
calculation shows that 



Condition ((TMjl <^ \A\ < 



(242) 



k' 



Condition ^fT^ <^ |A| < ■;^j^ . 



1 

2T2 



Since 2fc > k for all k = 1, 2, 3, ... , we see that if 
then condition 1(194(1 is satisfied but 1(195(1 is not. 
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Next, consider for M > and e > the example 
(244) V{0) = 






47r 2 ' ' 47r 2 



e< le*! <7r. 



On one hand, we have 



(245) v"{e) 



M(e-|0|)- \e\<e 

Ztt 



Me 
~2^ 



2 



e< |0| <7r, 



and so 

(246) Condition fBsl) ^ A/ < ^ 



On the other hand, direct calculation of the Fourier coefficients gives 
(247) Vk = 



2A/ , /fee , 



7Mt^ Me' Mttc^ 



k^O. 



487r 4 12 

Thus, we may estimate the sum in the condition p94|) as follows: 

^ ,,.., ^4Af ^fke\ ^;^ AM ,fke\ ^^^ Me^ 

fc=-oo fc=l ^ ^ A:=l ^ ^ fc=l 

since | sin(a;)| > \x\/2 for |a;| < tt/2. If the product Me^ is held fixed this lower bound can be made arbitrarily 
large (and in particular larger than one) simply by taking e sufficiently small due to the divergence of the 
harmonic series. We therefore see that if < c < Me^ < tt and if e is sufficiently small, then condition p95fl 
is satisfied but condition H194() is not. 
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